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Z usa m menfassu ng 



In dieser Arbeit wird die Theorie eines quantisierten Spin-2 Feldes behandelt. Dies 
geschieht im Rahmen der kausalen Storungstheorie nach Epstein und Glaser. Die Arbeit 
besteht aus zwei Teilen. Im ertsen Teil untersuchen wir die Eichstruktur eines mas- 
selosen selbstwechselwirkenden Spin-2 Feldes. Dabei nehmen wir einen reinen feldtheo- 
retischen Standpunkt ein, d.h. es werden keine geometrischen Aspekte der allgemeinen 
Relativitatstheorie vorausgesetzt. Aus dem Prinzip der Operatoreichinvarianz werden in 
erster Ordnung Storungstheorie notwendige und hinreichende Bedingungen abgeleitet, 
die eine eichinvariante Theorie eines Spin-2 Feldes erfullen muss. Dieses Prinzip besagt, 
dass die Eichvariation bzgl. der Eichladung Q der Selbstkopplung des Spin-2 Feldes eine 
Divergenz im Sinne der Vektoranalysis sein muss. Es zeigt sich, dass die allgemeinste 
trilineare Selbstkopplung des Spin-2 Feldes sich von der Einstein-Hilbert Kopplung nur 
um Co-Rander und Divergenzen unterscheidet. 

Im zweiten Teil dieser Arbeit (Kap.9) beschaftigen wir uns mit dem Langabstandsver- 
halten der Theorie eines Spin-2 Feldes welches an massive skalare Materie koppelt. 
Es wird der adiabatische Limes, bei dem die Abschaltung der Wechselwirkung im Un- 
endlichen entfernt wird, in Strahlungskorrekturen zur Zweiteilchenstreuung untersucht. 
Wir berechnen den Wirkungsquerschnitt fiir Graviton Bremsstrahlung in dem von einem 
der streuenden Teilchen ein Graviton niedriger Energie emittiert wird. Es zeigt sich, dass 
der Wirkungsquerschnitt im adiabatischen Limes logarithmisch divergiert. Desweiteren 
werden das Infrarotverhalten der Graviton Selbstenergie sowie die Selbstenergie der mas- 
siven skalaren Materie untersucht. Die Graviton Selbstenergie ist endlich im adiabati- 
schen Limes, wahrend bei der Selbstenergie der Materie ebenfalls eine logarithmische 
Divergenz im Wirkungsquerschnitt entsteht. 



Summary 



This work deals with the theory of a quantized spin-2 field in the framework of causal 
perturbation theory. It is divided into two parts. In the first part we analyze the gauge 
structure of a massless self-interacting quantum tensor field. We look at this theory 
from a pure field theoretical point of view without assuming any geometrical aspect 
from general relativity. To first order in the perturbation expansion of the S'-matrix we 
derive necessary and sufficient conditions for such a theory to be gauge invariant, by 
which we mean that the gauge variation of the self-coupling with respect to the gauge 
charge operator Q is a divergence in the sense of vector analysis. The most general 
trilinear self-coupling of the graviton field turns out to be the one derived from the 
Einstein-Hilbert action plus coboundaries and divergences. 

In the second part of this work (sect. 9) we consider a massive scalar field coupled 
to gravity. We are interested in the long range behaviour of this theory. Radiative 
corrections for two particle scattering are investigated in the adiabatic limit, where the 
cutoff of the interaction at infinity is removed. We compute the differential cross section 
for graviton bremsstrahlung in which one of the scattered particles emits a graviton of 
low energy. It is shown that such processes are logarithmically divergent in the adiabatic 
limit. Furthermore we show that the differential cross section for two particle scattering 
with a graviton self-energy insertion is finite in the adiabatic limit while for matter 
self-energy it is logarithmically divergent, too. 
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Geschrieben steht: „im Anfang war das Wort!" 

Hier stock' ich schon! Wer hilft mir weiter fort? 

Ich kann das Wort so hoch unmoglich schdtzen, 

Ich mufi es anders iibersetzen, 

Wenn ich vom Geiste recht erleuchtet bin. 

Geschrieben steht: im Anfang war der Sinn. 

Bedenke wohl die erste Zeile, 

Dafi deine Feder sich nicht ubereile! 

1st es der Sinn, der alles wirkt und schafft? 

Es sollte stehn: im Anfang war die Kraft! 

Doch, auch indem ich dieses niederschreibe, 

Schon warnt mich was, dafi ich dabei nicht bleibe. 

Mir hilft der Geist! Auf einmal seh' ich Rat 

Und schreibe getrost: im Anfang war die Tat! 

Goethe, Faust I 



1 Introduction 

One of the most interesting unsolved problems in modern theoretical physics is 
to combine Einstein's general theory of relativity with the principles of quantum 
physics. It is experimentally tested to a high degree of accuracy that the classical 
description of matter breaks down on the atomic scale and that quantum correc- 
tions become important. For the theory of general relativity one expects that it 
is valid up to the scale of the Planck length, Ip which is the unique combination 
with the dimension of length of the fundamental constants of nature, the speed 
of light, c, Planck's constant, h and Newton's gravitational constant, G, namely 
Ip = {Gh/ c^Y^'^ ~ 10~^^ cm. Beyond this scale it is believed that effects of a full 
theory of quantum gravity come into play. Therefore one might ask the questions: 
How does gravity behave at microscopic scales? What does the spacetime look like at 
the order of the Planck length? To answer these questions a full theory of quantum 
gravity is needed. Although there are many different approaches to the problem, 
none of them succeeded in overcoming the fundamental problems arising. One of 
the main problems in quantum gravity is that we have to give a meaning to the 
quantized metric field. Since the metric defines the geometry of spacetime we are 
forced to explain what a quantized geometry should be. The postulate of micro- 
causality becomes meaningless because there is no longer a fixed background on 
which concepts like timelike, spacelike or lightlike objects can be defined. 

In this work we deal with the theory of a covariantly quantized tensor field 
in four-dimensional Minkowski spacetime. The idea to describe quantum gravity 
in this way goes back to the work of Rosenfeld, Fierz and Pauli [5^, 18, 34| in 



the thirties and it was further developed by Feynman and deWitt in the sixties 
1 17, ^, 1^, |ll|. In this approach the metric tensor field g^u is splitted into two 



parts, the constant background field which describes the causal structure of 
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Minkowski spacetime and a dynamical part /i^jy which will then be quantized in this 
spacetime. These quanta of the gravitational field, called gravitons in the following, 
have zero mass and spin 2. Altough pure quantum gravity is finite at one loop it was 
realized by Deser, van Nieuwenhuizen and Boulware [0, ^ as well as by t'Hooft and 
Veltman in the mid seventies that a theory of a spin-2 field coupled to matter 
is non-renormalizable. Later in the eighties Goroff and Sagnotti computed the 
divergences of pure quantum gravity at two loop level. After the discovery that such 
a spin-2 theory is non-renormalizable people lost their interest in it. Clearly such a 
theory has not much predictive power because of the proliferation of undetermined 
constants which appear in the perturbative calculations. Nevertheless it was shown 
recently by Grillo |25| in the calculation of quantum corrections to the Newtonian 
potential that normalization terms only affect the potential at the origin. Therefore 
the long range part of the potential remains untouched. This will also be the case 
in higher orders, see [37|. In addition to that we think that one can learn something 
about the gauge structure of quantum gravity by considering this approach. The 
theory of a quantized tensor field will be treated here in the framework of causal 
perturbation theory which has the advantage that no ultraviolet divergences occur 
due to the mathematically correct treatment of the distributions. 



2 Causal Construction of the S'-Matrix 

Causal perturbation theory goes back to ideas of Bogoliubov |^ and was carefully 
developed by Epstein and Glaser [16| in the seventies. Later Scharf |36| applied it 
successfully to QED. In this approach one considers the S-matrix as the central 
object. The idea is to fix the first order interaction and then to construct higher 
orders of the scattering-matrix S by induction using free fields only. Since the free 
quantum fields as basic objects are operator-valued distributions the scattering- 
matrix will be constructed as an operator-valued functional on some test function 
space. The most important ingredient for the inductive construction is the causality 
requirement, which roughly states that the scattering-matrix of a sum of two test 
functions factorizes if the supports of the test functions can be separated in time. 
As test function space one considers for convenience the Schwartz space of functions 
of rapid decrease, since this space is invariant under Fourier transformation. Then 
all expressions are well defined tempered distributions. In doing this we cut off the 
interaction at large but finite distances which is unphysical in most cases and has to 
be removed at the end. This is the so called adiabatic limit where we take the limit 
that the test function goes to a constant. In this way we investigate the long range 
behaviour of the theory. 

The starting point for the construction is the S-matrix given as a formal power 
series in the coupling constant of the theory 

S{g) = 1 + '^— / d'^xi . . .(fxnTn{xi, . . . ,Xn)9{xi) . . . g{xn) (2.1) 
-I IT-- J 

n=l 
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where is a test function from the Schwartz space 5(M^) of rapidly decreasing 
functions. The Tn{xi, . . . , x„) are time-ordered products of interaction Lagrangians. 
They are given in terms of Wick monomials of free field operators. These free fields 



are operator- valued distributions on Fock space [49, 42]. If an ordering of the argu- 



ments of T„ would be given, such that > X2 > ■ ■ ■ > then we could write 

Tnixi, ... ,Xn) = Ti{xi) . ..Ti{Xn) (2.2) 

where Ti defines the interaction to first order of perturbation theory. In general no 
such time-ordering is given and one has to construct T„ carefully to every order. 
The naive construction of the T„ with the help of the step-function is ill-defined 

Tnixu... ,Xn)= -2;°(2))---^(^°(n-l) -^°(n))^l(^<T(l))---^l(^^(n)) 

(2.3) 

because this exression involves pointwise products of distributions which are not 
defined a priori. The use of this time-ordering would led to the well-known UV- 
divergences. The Epstein-Glaser method gives an inductive construction of the time- 
ordered products which are well defined and free of UV-divergences. In this method 
the 5-matrix is constructed by causality. This means that if we take two test func- 
tions gi and g2 for which we can find a separation of the supports in time by a 
spacelike surface, i.e. supp((7i) > supp(52) then the 5-matrix should factorize ac- 
cording to 

Sigi+g2) = Sig,)S{g2). (2.4) 
This means for the time-ordered products 

Tn {xi, . . . , 3^n) — ^m(^l ? • • • ; 3^m)^n— m {Xm+1 > • • • j X^) (^•^) 

if {xi, . . . ,Xm} > {xm+1, ■ ■ ■ ,Xn}- This factorization property lies at the heart of 
the causal construction. Now we assume that all Tm for 2 < m < n — 1 are already 
constructed. One can show that they have the form 

T-m {xi, . . . , Xfji) — ^ ^ • Ofe (xi , . . . , XjTj) . (xi XjYi ) • • • 5 Xfn—1 X-m) (^•6) 
k 

where : Ok{xi, . . . product of normally ordered free field operators and 

tm\xi — Xm, ■ ■ ■ ,Xm-i — Xm) IS a translation invariant numerical distribution on 
M^. Then Epstein and Glaser have shown that is a well-defined operator-valued 
distribution. For the explicit construction of the n-th order we proceed as follows: 
First of all we construct auxiliary distributions R'^{xi, . . . ,Xn) and A'^{xi, . . . ,Xn) 
from all the known of lower order by carrying out all possible contractions be- 
tween field operators appearing in T^yi and using Wick's theorem to obtain normally 
ordered expressions according to 

{X) (2.7) 

P2 



3 



,X^)=J2TnAX)Tn-nAy,^n) (2-8) 

where Tm are appear in the expansion of S{g)~^ . The sum runs over all partitions of 
the set of points {xi, . . . , Xm} into two subsets X, Y so that X is not empty. Then 
we introduce the difference 

Dnixi, ... ,Xn) = R'ni^l, ■■■ ,Xn) - A'^{xi, ... , Xn) (2.9) 

which has the following form 

Dnixi,. .. ,x„) = ^ : Ok{xi, ... ,Xn) : dl^\xi - x„, . . . - x„) (2.10) 

k 

where : Ok{xi, . . . normally ordered product of free field operators and 

dih\xi — Xn, . . . ,Xn-i — Xn) is a numerical distribution that is build out of prod- 
ucts of positive and negative frequency parts of the Pauli-Jordan distribution. Such 

(k) 

products are well-defined. Due to the translation invariance, dn depends only on 
n — 1 independent coordinates. The most important property of the distribution Dn 
is causality, i.e. 

supp(i:'„(xi, . . . ,xn)) ^ r+_i(x„) ur;;_i(x„) (2.11) 

where 

r^_i(xn) = {(xi,... ,x„) eM4"|xj G (x„ + y±),Vj = 1,... ,n-l}. (2.12) 

The support properties are entirely encoded in the numerical part of Dn. The next 
step in the construction of T„ is to split the distribution Dn into retarded and 
advanced parts, Rn and An, respectively such that 

supp(i?„(xi, . . . ,Xn)) ^r+_i(x„) 

supp(A„(xi, . . . ,xn)) ^r;;_i(x„). 

(k) 

To achieve this we have to split the numerical part dn . The critical point for this 
operation is the total diagonal, i.e. 

r+_i(x„) n r;_i(x„) = a„ = {xi = . . . = x„}. (2.13) 

(k) 

Due to the fact that d„ is translation invariant we can shift the critical point 

(k) 

to the origin. The behaviour of the distribution dn at the origin is measured by 
the singular order. The singular order can be defined for an arbitrary distribution 

Definition 1. One says that the distribution t has scaling degree s at x = 0, if 

s = inf{s' G M|A' t{Xx) ^ in the sense of distributions}. 

In spacetime dimension d = 4 the singular order is then related to the scaling degree 
by uj := [s] — 4, where [s] is the greatest integer less or equal to s. 
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In order to avoid UV-divergences one has to do the sphtting operation with the 
correct singular order. It turns out that if the singular order uj < then the splitting 
is trivial and the multiplication with the discontinuous step function can be done 
without the appearence of UV-divergences. On the other hand, if u; > then the 
splitting operation is non-trivial and moreover non-unique. In this case we obtain 
for the retarded part 

^ ("^1 -^m • • • ) ^n—l ^n) ^ ^ (''^1 •^m ■ ■ ■ : ^n—l ^n) 

UJ 

+ Yl Ca,kD%^*'^''~^^^ (Xi - X„, . . . , Xn-l - Xn) 
\a\=0 

(2.14) 

(k) (k) 

where w > is the singular order of the distribution dn and rn is a special 
splitting solution. The splitting solution is therefore only determined up to a sum 
of normalization terms with finite coefficients Ca^k- They are not determined by the 
splitting procedure itself and there must be imposed further physical conditions like 
Lorentz covariance, gauge invariance, etc. to restrict them. The retarded part can be 
determined in momentum space by a dispersion integral ||3^ . Finally r„, including 
it's normalization terms is given by 

Tn{xi 1 • • • 1 — Rn (^^l ; • • • > Xn) (xi , . . . , Xn) 

— ^ ^ • Of^ (xi , . . . , Xn) ■ ^ (-^1 Xn j • • • j Xn—1 Xn) ■ ^ ^ 
k 

In this way we obtain a well-defined time-ordered n-point function which is renormal- 
ized without introducing any regularization presrciption. With this T„ it is possible 
to compute UV-finite physical quantities, like amplitudes or cross-sections. 



3 Gauge Invariance in Classical Linearized Gravity 

The general theory of relativity can be derived from the Einstein-Hilbert Lagrangian 

Leh = -^V^R (3.1) 



where R = g^'^R^u is the Ricci scalar, g is the determinant of the matrix {g^'^) and 
= 327rG {G is Newton's gravitational constant) |45, 48]. It is convenient to work 
with Goldberg variables |19| 



9' 



-99 



(3.2) 



As was already mentioned in the introduction we will in this section consider the 
linearized theory of gravity, therefore we expand the inverse metric g'^'^ in an asymp- 
totically flat geometry 



(3.3) 
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Here rj'^'^ = diag(+l, — 1, — 1, — 1) is the metric of Minkowski spacetime. All tensor 
indices will be raised and lowered with r/^'^. The quantity h^'^ is a symmetric sec- 
ond rank tensor field, which describes gravitons after quantization. Formally (3.1) 
becomes an infinite power series in k 

oo 

Leh = ^'LeI (3-4) 

j=0 

The lowest order term is quadratic in h^'^{x) and defines the free asymptotic 
fields. The linearized Euler-Lagrange equations of motion for h^^{x) are 



□/i^^(x) - -r]^'''Uh{x) - h^'f^^ix) - Kf^^^^x) = (3.5) 



where h(x) = h^{x). This equation is invariant w.r.t. gauge transformations of the 
form 

h^"" — > h^'^ + u^'''' + n'^'^ - 7]'"'uPp (3.6) 

where is a vector field which satisfies the wave equation 

□n^(x) = 0. (3.7) 

These gauge transformation correspond to the general covariance in it's linearized 
form of the metric tensor g^y[x). The corresponding gauge condition, compatible 
with (|3.6|) is the Hilbert-gauge 



h^^ = 0. (3.8) 
Then the dynamical equation for the graviton field h^^'^ reduces to the wave equation 

Uh^''{x) = 0. (3.9) 

The first order term L^^}j gives the trilinear self-coupling of the gravitons 

Lfn = \ h'" {hfhf - \ h,^K + 2 K^hl"^ + /i>.<" - 2 K^K^) . (3.10) 

There exists many alternative derivations of this result ( ^.10 ), starting from 
massless tensor fields and requiring consistency with gauge invariance in some sense 
1 44, 32, 33]. In classical theory the work closest to our non-geometrical point of 
view is the one of Ogievetsky and Polubarinov [^]. These authors analyze spin- 
2 theories by working with a generalized Hilbert-gauge condition to exclude the 
spin one part from the outset. They impose an invariance under infinitesimal gauge 
transformations of the form 

J/i/^^ = d'^u" + d'^u'' + rj'^'dau'^ (3.11) 
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and get Einstein's theory at the end. Instead W yss [50| ] considers the couphng to 
matter. Then the self-couphng of the tensor-field (|3.10| ) is necessary for consistency. 
Wald [^] derives a divergence identity which is equivalent to an infinitesimal gauge 
invariance of the theory. Einstein's theory is the only non-trivial solution of this 
identity. In quantum theory the problem was studied by Boulware and Deser 
These authors require Ward identities associated with the graviton propagator to 
implement gauge invariance. All authors get Einstein's theory as the unique classical 
limit if the theory is purely spin two without a spin one admixture. 



4 Perturbative Quantum Gauge Invariance 

In this work we will study the problem without any reference to the metric tensor by 
means of perturbative quantum gauge invariance. This method which was worked 
out for spin-1 non-abelian gauge theories (massless [P and massive []T5| , ^) in last 
years proceeds as follows: First one defines infinitesimal gauge variations on free 
fields. In the case of tensor fields it looks like ( 3.1l| ) where u^{x), instead of being 



an arbitrary function, is now a Fermi field which satisfies the wave equation. u^{x) 
may be regarded as a free Fadeev-Popov ghost field. These ghost fields play a very 
important role in connection with gauge invariance. We write down the most general 
trilinear coupling Ti between the graviton and ghost fields which is compatible with 
Lorentz covariance, power counting and certain basic properties (like zero ghost 
number). Next we impose first order gauge invariance which strongly restricts the 
form of Ti. Among the possible solutions we recover Einstein's theory L^eIi- The 
general solution can be written as a linear combination of L^^]j and divergences 
as well as coboundaries. In the perturbative construction of the S'-matrix we next 
have to calculate the time-ordered product T{Ti{x)Ti{y)} = T2{x,y) by means of 
causality ||l^, Then Schorn has shown that second order gauge invariance 
gives further restrictions, in particular, in the case of gravity it requires quartic 
normalization terms of the form L^^h iii the above expansion ( |3.4D . In this way the so- 
called proliferation of couplings can be overcome by perturbative gauge invariance. 

Our fundamental free asymptotic fields are a symmetric tensor field of rank 
two h^^"{x) and ghost and anti-ghost fields u'^{x) and u'^{x). We consider these 
fields in the background of Minkowski spacetime. A symmetrical tensor field has 
ten degrees of freedom, which are more than the five independent components of a 
spin-2 field. The additional degrees of freedom can be eliminated by imposing two 
further conditions p^, namely 



h^"'{x),^ = and /i%(x) = 0. (4.1) 
They are disregarded in the construction of the gauge theory and must be considered 



later in the characterization of physical states |22, 23|. 

Our tensor field h^'^{x) will be quantized as a massless field |2^, 25] as follows 

[/i"^(x), h^"'{y)] = -^rf^^"'Do{x - y) (4.2) 
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where Dq{x — y) is the massless Pauh- Jordan distribution and the tensor b"'^'^^ is 
constructed from the Minkowski metric r/^'^ in the following way 

l^aPf^u ^ 1 (j^"M^/3^ + ^^'^^f^f^ _ jfl^ri^^^ . (4.3) 
We can write down the Fourier-representation of the field W'^ix). It is given by 

= (27r)-3/2 I (a°^(fc) exp(-ifcx) + a"^(fe)t exp{+ikx)) . (4.4) 



Here u!{k) = \k\ and a°'^{k), a'^^^k)^ are annihilation and creation operators on a 
bosonic Fock-space. From (4^) we find that they have the following commutation 
relations 

[a"^(fc),a^^(fe')1 =6"^^^5(3)(^-^'). (4.5) 
In analogy to spin-1 theories one introduces a gauge charge operator by 

Q = [ h'^^{x),pdQu'^d^x. (4.6) 



For the construction of the physical subspace and in order to prove the unitarity of 
the S'-matrix we want to have a nilpotent operator Q. Therefore we have to quantize 
the ghost fields with anticommutators 

{u^^{x),u-{y)]=ir^^'''D^{x-y) (4.7) 

and all other anti-commutators vanishing. All asymptotic fields fulfil the wave equa- 
tion 

□ /i'^^(x) = 

□ ^i"(x) = (4.8) 
Uu'^{x) = 

The gauge charge Q (|4.6| ) defines a gauge variation by 

dQF:=QF-(-l)"«(^)FQ (4.9) 

where Ug is the ghostnumber. This is the number of ghost fields minus the number 
of anti-ghost fields in the Wick monomial F. The operator dq obeys the Leibniz rule 

dQ{AB) = {dQA)B + (-l)"«(^)^dQB (4.10) 

where A and B are arbitrary operators. We obtain the following gauge variations of 
the fundamental fields: 

rfQ^^' = -^K + <M-<'^<) (4-11) 
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dQh= iut^ (4.12) 
dQui'= ihif^ (4.13) 
dQU^'= (4.14) 

From ( 4.11J ) we immediately see 

dQ/if/ = 0. (4.15) 

The result ( |4TlD agrees with the infinitesimal gauge transformations of the Goldberg 
variables, so that our quantization (|4.2| ) and choice of Q corresponds to the classical 
framework described in section 3. The asymptotic fields will be used to construct 
the time-ordered products T„ in the adiabatically switched S-matrix 

S{g) = 1 + ^^— jd^xi...d'^XnTn{xi,...,Xn)g{xi)...g{xn) (4.16) 

n=l 

where g G 5(R^) is a test function. The time-ordered products T„ are operator- 
valued distributions and they can be expressed by normally ordered products of free 
fields, see section 2. It is very important that gauge invariance of the S'-matrix can 
be directly formulated in terms of the T„ [37|. First order gauge invariance means 
that dqTi is a divergence in the sense of vector analysis, i.e. 

dQT^{x) = id^T^^^{x). (4.17) 

The definition of the n-th order gauge invariance then reads 

^ Q 

dqTn = [Q,Tn] = i^T-p:T^/i{xi,. . . ,xi,... ,x„). (4.18) 
1=1 

Here T^^^ is the time ordered product with a gauge variated vertex Ty^{xi) at 
position xi and ordinary vertices Ti at the other arguments. 



5 Structure of the Interaction 

Here we introduce the self-couplings of the quantum tensor field h^^ (x) . The simplest 
expression leading to a self-interacting spin-2 field theory is a trilinear coupling of 
the quantum fields h^'^{x) and h{x) = h^^{x). We require Lorentz invariance and in 
addition to that two derivatives acting on the fields. This is for the following reasons: 
First of all, by inspection of all trilinear self-interaction terms without derivatives, it 
is easily seen that such a theory cannot be gauge invariant to first order of perturba- 
tion theory. Therefore an interaction without derivatives can be ruled out. Secondly, 
it is impossible to form a Lorentz-scalar from three rank-2 tensor fields with only one 
derivative. Last but not least the corresponding trilinear expression in the expansion 
of the Einstein-Hilbert Lagrangian contains two derivatives as well. Therefore we're 
able to reproduce the results from classical general relativity. 
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(5.1) 



In the following all fields are free fields obeying the free field equations of motion. 
All products of two or more fields at the same spacetime point x are viewed as normal 
products. Then the general ansatz for a combination of three field operators contains 
12 terms[|: 

ri^(x) := ai : /i^'^(x),^/i(x),^/i(x) : +02 : h^'^ {x)h{x),^h{x),u : 

+ 03 : : +04 : /i"'^(x)«/i^^(3;)/i(x),^ : 

+ 05 : /i"'3(x)/i^'^(x)«/i(x),^ : +06 : /i"'3(x),^/i^'^(x),a/i(x) : 

+ 09 : h'"'{x),ah''''{x),,3h^'^{x) : +aio : h^"" {x)„h''''{x)h^'^ {x),p : 
+ aii : h^"'{x)h'"'{x)^h^'^{x),p : +ai2 : h^'^ {x)h''''{x),ph^'^ {x)^ : 

Here we have omited all terms which are divergences. These are terms with a con- 
traction on the two derivatives, e.g. 

h{x),ah{x),ah{x) = 1/2 da [h{x),ah{x)h{x)Y 

Furthermore all terms with two derivatives acting on the same field can be trans- 
formed into a divergence plus a term already contained in ( |5.1D . These terms would 
modify our ansatz only in a redefinition of some parameters and can be omited 
without losing generality. 

As in the cases of Yang-Mills theory |l^, ^] and Einstein gravity |39| we expect 
to get a gauge invariant first order coupling only if we couple the tensor field h^'^ also 
to ghost and anti-ghost fields. The most general expression with zero ghost-number 
is 



T-ix) 



bi : uP{x),^u''{x),ph'"'{x) : +62 : u^(x),^ti^(x)/i'"^(x),p : 
+ 63 : uP{x)u^'{x),^h^"'{x),p : +64 : uf{x),pU^'ix),^hf"'{x) 
+ 65 : uP{x),pU^'{x)h^''{x)„ : +^6 : uP{x)u^{x),ph^"' {x)„ 
+ bj : uP{x),pU^{x),ph{x) : +^8 : uP{x),^u^{x)h{x),p : 
+ 69 : uP{x)u^'-{x),f,h{x),p : +610 : (x) ,pU^' (x) ,f,h{x) : 



+ 611 : uP{x),pu''{x)h{x) 
+ hi^:uP{x),^u^'{x),,hP%x) 
+ bi^:uP{x)u^'{x),phP''{x),, 
+ biT.uP{x),,u''ix)hP^x),p 
+ 619 :u^(x),^n'^(x),p/i'"^(x) 
+ b2i:uP{x)uf'{x)„hP''{x),p 



,p : +612 : uP{x)uP{x),ph{x),p, : 
+ bu:uP{x),f,u''{x)hP''{x),^ 
+ bie:uP{x)„uP{x),^hP''{x) 
+ bis:uPix)uPix)„hP''ix),^ 
+ b2o:uPix),^uPix)hP''ix),p 



(5.2) 



^We use the following convention regarding the indices. All vector and tensor indices are written 
as superscript, whereas all partial derivatives are written as subscript in the abbreviated form with 
a prime in front of the index, i.e.: A{x),i, = dA(x)/dx". All indices will be raised and lowered by 
the Minkowski metric ly^j^ and will be properly contracted like A'^B'^ ~ rji_ivA'^B^ . 
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We will suppress all arguments of the field operators as well as the double dots of 
normal ordering in subsequent expressions. The complete first order coupling is then 
given by: 

ri(x) :=r/^(x)+r«(x). (5.3) 

In the following analysis we want to study in which way the parameters of the theory 
ai, . . . , ai2 and 61, . . . ,621 will be restricted due to first order gauge invariance. 

6 Classification of Divergences 

In the previous section we have defined the trilincar coupling of the graviton field 
Ti as well as the coupling to ghost and anti-ghost fields T". In this section we try to 
write the gauge variation dgTi as a divergence d^Ty^. We proceed in a systematic 
way: Because of the great variety of different terms in dqTi it is most convenient to 
use a separate ansatz for T^^^ . Since the operator dq applied to our Ti increases the 
ghost number of the result by one we have to make an ansatz with ng{T^^-^) = 1. 
Furthermore the application of dg increases the number of partial derivatives by 
one. Due to this, every term in must have two derivatives. The terms appearing 
in this ansatz can be classified according to their field content. In the so called 
graviton sector we have Tj^^ ~ uhh and the ghost sector contains Tj^^ ~ uuu. 
Inside each sector there is a further classification w.r.t. the tensor indices: There are 
seven different Lorentz types in the graviton sector, namely 

type A : uVhf^ 
type B : u^hh 
typeC : u"/i"^/i 
type LI : u'^h'^^h^P 
typcL : u'^h'^'^hP'' 
typeF: ti^/i^'^/i^'^ 
typeG : vl'h'^h. 

In the ghost sector we've three different Lorentz types, namely 

type if : 
type J : u°'u'^uP 
typeK : u^u^u". 

In the following subsections we explain in detail the way in which the divergences 
for these different Lorentz types can be found. 



11 



6.1 Graviton Sector 

Divergences of type A have the structure dadudi,\u°'hP^h'"^ , where we have to dis- 
tribute the three derivatives in all possible ways among the fields. Taking into ac- 
count that all fields obey the wave equation we find the so called basis elements from 
which all divergences of type A can be constructed: 



For the construction of the divergences we have two different partial derivatives with 
indices a and v which can be taken out, i.e. we can form 



where the remaining derivatives inside the bracket must be distributed in all possible 
ways among the fields. Then we find the following divergences of the form div^ 



d,{ulhPJhP'') =et + ei 
d,{uXhP'')=ei + ei 
d,{u"hP^„hP'')=ei + ei 
d,{u'^hf:^K")=ei + et 



In this way we've found the set of different divergences of Lorentz type A. Collecting 
all terms in the brackets we can make the following ansatz for Tt^jf': 



rf;f = ci uXh'" + C2 ui'hP^hP^ + C3 u'^,^hP''hP'' + C4 
+ C5 uZhP^hP'' + C6 u'^hP^hP^ + cr uXhP'' 




divi 
div2 



d^{dad:,\u'^hP"hP") 



d^{uZK"hP")=e^ + e^+ei 
da{u''hP^hP,^)=ei + 2ei. 



Divergences of the form div. 



'2 are 
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where the constants ci , . . . ,07 are for the moment free constants to be determined 
by gauge invar iance. 

Divergences of type B have the structure dad,yd,^\u"hh. The basis elements of 
this type are 



e? 


— v" h h 


ef 


— ?/"/? h 


ef 


— y"h h 


ef 


— v"h h 


ei 





Again we can form two different types of divergences corresponding to the two partial 
derivatives with indices a and v. 

divf = da{di,di,\u°'hh) 
div2 = di,[dadi,\u°'hh) . 

We find the following divergences of the form divf 

d4uZh,,h) =ef + ef + ef 
5aK^,A.) =ef + 2ef. 

For divergences of the form div2 we find 

du{ul,,hh) = 2ef 

=ef + ef 

For T[yj^ we can therefore make the ansatz 

T'^jf = C8 </i«/i + eg ui^h^h^ + Clo <,^/i/i + cii u°'h^,^h 

+ C12 uyi,Jl + Ci3 U°'hjl,^ + Ci4 u^i^h. 

Divergences of type C have the structure dfjidi,di,\u°'h°'^h. The basis elements of 
this type are 

pC _ a uafii 



en = u' 
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'^tv"'^fl "''1/ 










4 





We find two types of divergences 

div^ = d^{di,d^\u''h''^'h). 

We find the following divergences of the form divi 

d^{ulhrh)=e^ + (^,+e^ 

d^{u-hrK)=e'i + e'i + e^. 
For divergences of the form div2 we find 





= e? 


^4 


d,{u%hrh) 


= ef 






= e? 


+ ef 






^4 


d,.{uZh^^h,^) 


= e^ 


+ 4 




= e? 


+ 4 




= e^ 


+ 4 


d,[u^hZ^h,^) 


= e? 


+ 4 




= 4 


+ 4 



For T^'^ we can make the ansatz 

rf/'f = ci5 uZhZ^'h + ci6 uZ,h''^'K + ci7 + cis 

+ Ci9 n^/l^';/! + C20 U'^h'^'h,,,^ + C21 + C22 (6-3) 

+ C23 + </ir^ + C25 + C26 

Divergences of type D have the structure dpdydi,\u'^h"-'^h'^P . The basis elements 
of this type are 
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e? 




ef 


= u,i,h,p,i,h ^ 


ef 




ef 


— v°'h°"^h'^P 


ef 




ef 




ef 





We find two types of divergences 

dwf = 5^(ap5^|n°/i"^/i''''). 

We find tlie following divergences of type divf' 

dp{uZhZ-h''P)=e? + ef + e^ 
5p(0°-/tf/) =ef + ef + ef 
5p(u"/tr^f/) =ef + ef + ef. 

Divergences of type cZwf are 





= ef 


+ ef 




= ef 


+ ef 




= ef 


+ ef 




= ef 


+ ef 




= ef 


+ ef 


d,{u''h%lh''c) 


= ef 


+ ef 


d,{u''h'^;hlP) 


= ef 


+ ef 




= ef 


+ ef 




D 1 D 
= 65 +69 



For Ty-^ we can make the ansatz 

Tf/'f = C27 U°/t«'^/t-'^ + C28 + C29 U^/l^/lf/ + C30 U%,^h'^'' h'^P 

+ C31 u'^hUh'^P + C32 u'^h'^^h'^P^ + C33 ^^/i^'^/i'^'' + C34 

+ C35 tX^/l^/lf/ + C36 + C37 + C38 ^^^/ij'/ifp". 
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Divergences of type E have the structure dpdadi,\u°'h°'^hP'^ . The basis elements 

are 



Although there arc three different partial derivatives for type E we have only two 
different divergences because of the symmetry of the tensor field h'"^ 



divf = dp{dad^\u"h"''h^'') 
divf = d^{dpda\u"h'"'h''''). 



We find the following divergences of the form div 



dp{u^hZlhP^)=ei + ef, + ef, 




u,ph ^CT'i/ 
'^ti/h,p 
^ h,p^,i/h^ 

u h,y h^p-fj 



=ef + ef + ef 



16 



Divergences of the form div2 are 

d,{u?^hZ''hf"')=ef + ei + ef 
d,{u%h''-hP^)=ef + el + el 

d4n'^h'^''h^p=eg + ef, + ef,. 
For we can make the ansatz 

Tf/'f = C39 uZ,ph'^Ph^''' + C40 U'^h%%hl"' + C41 u'^h'^Ph'^.^ + C42 t^^/l^^'/t'^'^ 
+ C43 + C44 + C45 tX^/t^/l'^'^ + C46 

+ C47 t^^/l^/l^r + C48 <p^h'^^'h''P + C49 n^/l^^^/i'"^ + C50 n°/t"'^/if/^ 
+ C51 tX^/l^'^/t'"^ + C52 + C53 W^K^hP^. 

Divergences of type F are of the form dpdadnlu'^hP^hP''^ . The basis elements are 



ef 






ef 












ef 






ef 




^jph^,i/h^ 






^'ph^^^KT-v 


ef 




^'p^^^n/ 








ef 




"^tl/^^pKT^^ 


^10 




"^tl/^^p^^V^ 














^13 




^ h^ptu^w 


F 

ei4 






ef5 
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Then we have two different types divergences, namely 

divf = dp{d^du\u''hP''h'"') 
div^ = d^{dpd^\u''hP''h''''). 

We find the following divergences of the form div[ 

5,«^«^/i^'^)=ef + ef + ef 
a,«mr)=ef + ef + ef 

5,(u-^«'^)=ef + ef3 + ef5 

Divergences of the form div2 are 

a.Kp^/i^^/i'^'^) =ef + 2e^ 

=efo + 2ef4 

For Ty^ we can make the ansatz 

Tf/'f = C54 u''„,X'h'''' + C55 + C56 n^/l^'^/ir^Tp + C57 ^f^f/Zl'^'^ 

+ C62 n^'/lf.'^/ir." + C63 + C64 u'' h^p'^h'"' + C65 ^x'^/l^T/ir/ 

+ C66 ^/iSr^"" + C67 ^/i'^'C" + C68 U'^h^^'h-'J. 

Divergences of type G are of the form dpd,jdi,\u^h'"^h. The basis elements 

&\ = u^p,(T'uh^'^h 

64 ^ip^u^KT ^ 
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^5 


= 






= 


^^p^rr'i/^ 




= 




^8 


= 


"'ip"'^ "''V 


69 


= 


"''p"'v/ "'rr 


^10 


= 




p° 
^11 


= 










^13 










^ ^^pXT^'V 






^ ^^pn/^KT 


^16 




u h^p h,^,i, 








ef8 







We have two types of divergences, namely 

divf = dp{dadu\vyh'"'h) 
div^ = d,,{dpda\vrhP''h). 

We find the following divergences of type divf 



dpiui^^h^^h) 


= ef 


+ + ef^ 


dp{uihi^h) 


= ef 


+ + 


dp{uihf"'K) 


= ef 


+ + 






+ eg + efs 


dp{u-h^K) 


— '^8 


+ + efg 


dpKhP^^h) 


— 64 


+ efg + 


dp{u'^hp^h„.^) 


- P^ 

— '^7 


+ efe + efs 


dp{u'^hi-K^) 


- P^ 


+ efs + efr 


dp{uihp^h^) 


- p"^ 


+ + efi 



Divergences of the form div2 are 

=pf + ef + ef 
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d,{u''hP-h,p^)=ef, + ef, + eg. 
For we can make the ansatz 

r{)'f = C69 + C70 n^'/i^';,^ + C71 + C72 

+ C73 uP,h'"'h,p + C74 ■u^'/i^^/i,^ + C75 nf^/i(^''^ + C76 uPph>^h^ 
+ C77 + C78 U^^hP'^h + C79 n'^/lf/^/t + C80 

+ C81 nf^/i^'^/l + C82 nf^/t'"^/t^ + C83 u^^hi^h,^. 

This completes the discussion of the divergences in the graviton sector. 
6.2 Ghost Sector 

Now we come to the ghost sector where we've three different Lorentz types. Diver- 
gences of type H are of the form dadvdu\u°^u'^u°^ . The basis elements are 



^1 






p" 






p" 






p^ 






p" 






p" 




<« 


p" 




<« 


p" 






^9 







We have two types of divergences corresponding to the two different partial deriva- 
tives 

div^ = d^{dad^\vf'u''u''). 

We find the following divergences of the form div^ 

=ef + ef + ef 
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Divergences of the form div2 are 
Then we can make the ansatz for Tj^jf 

T^f = C84 + C85 <U"< + C86 U^^^ + Cgr ^U^'^ 

+ C88 + C89 + C90 + C91 (6-8) 

+ C92 u^uf^u^ + C93 u^,ij_u°'u"' + C94 u''uZ,fj,u"' + C95 u'^u'^u'^,^. 

Divergences of type J are of the form dadpdv\u'^uPu^ . The fact that the ghost 
fields anticommute reduces the number of independent basis elements considerably, 
e.g terms of the form u'^,^.^uPu'^ , u'^u^puf^ and u'^u^uufp must vanish. Then we have 
the following basis elements 









e{ 




■~Ct p V 


4 






4 






4 






4 






4 




"''p"''a"'i' 


4 




"',p"',i/"'ta 


69 




"'VX'P'l/"' 


4o 










yOt p V 


^12 




yOt p V 

"'tptu"''a- 



We have two types of divergences, namely 

divf = da{dpd^\u''uPu'') 
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We find the following divergences of type divf 

da{u?pKu") = ei + ei + ei 

da(u°'u'!p,uu'') = 64 + eg + ei2. 
Divergences of type div2 are 

dpiu'^^u^u'C) = ei - ei 
dp{u''u^^,X) = 4 + 4 + 4i 

d^{uZu^u'') = ei - ei - ei 
dp{u'^uPu^^,^) =ei-eio-ei 

dp{uZuPu':,) = -4-4-4. 
For Ty"^ we can make the ansatz 

Tf/'/ = C96 n^n^n" + C97 uluf^uP + cgg + C99 nf^,^ii"n'' 

+ cioo uPu^X + cioi uPu"u^:„p + C102 <ii>'' + C103 nf„n"tx(; (6.9) 

+ C104 uPu^yUtp + C105 UPpU^yU^ + Ci06 u'^^U^'u'i^ + Ci07 U^U^U^^. 

Divergences of type K are of the form dad,ydi,\u"u'^u'^ . The basis elements are 



ef 


— ?7" ?y'^i/'^ 


ef 




ef 


""11/ "'i^ 


ef 





Then we have two types of divergences 

div^ = dc{dydy\u'^u''u'') 
div^ = d^{do,d^\u'^u''u''). 

We find the following divergence of type divf 

=ef + ef -ef. 
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K 



The following divergences are of type div2 



e 



„K 
1 



9,(^2"<„n'^) =ef + ef 
5,(7i"<nX) =ef + ef. 

We can make the ansatz for T^jf 

^l/f = ^108 u'^a^'itU" + C109 Ka'uU^u'' + CnO U^^U^ u" 

' (6.10) 

+ Cm u^n^n;^ + cii2 u'^^u'^u'' . 

This completes the discussion of the ghost sector. Collecting all divergences we 
obtain the total divergence as the sum of these 10 different Lorentz types 

W)i = ^M E ^i/v (6-11) 

i&{A,...,K} 

The parameters ci, . . . , cii2 G C are for the moment free constants, to be determined 
by gauge invariance. This expression^ for Ty-^ contains all possible combinations of 
fields appearing after gauge variation of Ti . Without losing generality one can now 
eliminate a few terms in the types A, . . . ,D,H and K^. For that purpose we consider 
a new Q-vertex T^^^{x) for which the following relation holds 

r{)^(x) = ff/,(:r) + S'^(x) (6.12) 

where has the special form B^{x) = df,A^^{x) and A^^{x) is an anti-symmetrical 
tensor of rank 2. Then we have 

9^r()^(x) = a^f{)^(x), (6.13) 

because partial derivatives are commuting. Let us now construct such a tensor A^^^. 
We consider the type-j4 term C3 u'^^^h^'^ h^^ . This can be written as 

C3 ul^^hP'^hP'^ = C3 [d^ {u-^^hP^hP^) - 2 uT^hP^hP''] . (6.14) 

In an analogous way and using the wave equation we can write 

= C3 uP^^hP-hP- = C3 [d^ {uP^hP'^hP'^) - 2 uXh"''] ■ (6.15) 

Now we add — C3 u^,ahP'^hP'^ to Ty-^ and obtain 

^Vi = % + C3 5. (n; - O hP'^hP^ . (6.16) 

■^T^'y-^ is called Q-vertex in the sequel because it is obtained from the usual vertex Ti if one 
replaces a quantum field with the gauge variation of that field. 
^This relies on an idea of M. Diitsch, see [p^ . 
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The expression in brackets is anti-symmetric in z/, /i and we get if we replace the 
constants ci with ci + 2 C3 and C7 with C7 — 2 C3 in T^^- In this way we can ehminate 
the monomials with constants C3, C4 in type A, cio, cn in type B, cig, cig, C20 in type 
C) C30, C31, C32 in type 093,094, C95 in type H and C109 in type K. Then we obtain 
a smaller Q- vertex Tj^^ from if we replace 

Q, i e {1, 2, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 21, 22, 23, 24, 25, 26, 27, 28, 29, 
33, 34, 35, 36, 37, 84, 85, 86, 87, 88, 89, 90, 91, 92, 108, 111, 112} 

by 

ci = ci + 2 C3 + C4, C2 = C2 + C4, C5 = C5 - C4, ce = ce - C4, 

C7 = C7-2C3, C8 = C8 + 2cio + Cii, Cg = Cg + Cu, C12 = Cu - Cn, 
Cl3 = Ci3 - Cii, Ci4 = Ci4 - 2cio, C15 = C15 + C18 + C19, C16 = C16 + C18 + C20, 
Cl7 = Cl7 + Cl9 + C20, C21 = C21 — Cl9, C22 = C22 — C20, C23 = C23 — C20, 
C24 = C24 - C18, C25 = C25 - C18, C26 = C26 - C19, C27 = C27 + C30 + C31, 
C28 = C28 + C30 + C32, C29 = C29 + C31 + C32, C33 = C33 - C31, C34 = C34 - C32, 
C35 = C35 - C32, C36 = C36 - C30, C37 = C37 - C30, C38 = C38 - C31, 

C84 = C84 + C93 + C94, C85 = C85 + C93 + C95, C86 = C86 + C94 + C95, C87 = C87 — C94, 
C88 = C88 — C95, C89 = C89 — C95, C90 = C90 — C93, C91 = C91 — C93, 

C92 = C92 — C94, C108 = C108 + C109, Cm = Cm + cio9, cu2 = C112 — ciog. 

In the following we will always use this new Q-vertex Ty^. After elimination of these 
redundant terms in the types A, . . . ,D,H and K one can express the corresponding 
terms of (IqTi in an unique way as a divergence in the sense of vector analysis. This 
is done in appendix A. For the types E, F, G and J the situation is different. Here we 
have only monomials without derivatives acting with respect to x^. For these types 
we obtain a Q-vertex which contains free constants. But it is shown in appendix B 
that this indeterminancy can be put into a form which drops out when we build the 
divergence. 

7 Relations from First Order Gauge Invariance 

Now we will work out the consequences of first order gauge invariance. This means 
that we compute the gauge variation of the first order coupling Ti and require that 
the result can be written in the form d^Ty-^, where is given by (|]l|) with the 
redundant terms eliminated according to the discussion above. 

7.1 Type E Divergences 

In this subsection we consider the type E divergences explicitly. From the comparison 
of these divergences with dQTilxypeE we will get linear relations among the coupling 
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parameters of Ti. We require the following equation to be satisfied 

dgTilTypeE = d^Tyf. 
Calculating the right side of this equation we get 



(7.1) 



The new constants are defined as follows 



C?29 
^32 
d35 

dss 
dii 
diA 



C42 + C48 , 
C43 + C48 , 
C44 + C49 , 
C45 + C49 , 
C46 + C50 , 
C39 + C48, 



C^30 

^36 
^9 

^45 



C40 + C42 + C51, 
C43 + C44 + C51, 
C41 + C44 + C53, 
C45 + C46 + C53 , 
C40 + C47 + C53, 
C40 + C49 , ^46 



From equation ([7.ip we see that 



C?31 
C^34 
C^37 
^40 

^43 



C42 + C47 + C52 
C41 + C43 + C52 
C39 + C45 + C51 
C39 + C46 + C52 
C47 + C50 



C4I + C50. 



(7.2) 



(7.3) 



C?29 
^32 
^35 

^38 

dil 
du 



2«io, 



0, ^36 



^30 = 
C^33 = - 

= -ibis, 



-ibi9, dsi 
■iai2, ^34 = 



-i{aii + 614 + 620) 



-ib 
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d 



37 



"M 2a8 + 2 "9 



0, ^39 = 

-i6i6, di2 = -ib2i, ^43 
0, ^45 = 0, diQ = 0. 



07 + ^17), d^o = -i[a7 + -aio 



(7.4) 



Finally we arrive at the divergence form if we invert the system (^^). This is done 



in appendix B. Let Me G Mat{18 x 15, Z) be the coefficient matrix of (7.3). Then 
we can write this system of equations as 



(7.5) 



where S C^^ and G C^^ are the column vectors with components (C39, . . . , C53) 
and (^29 5 • • • ,'^46) respectively. Now we observe two things: 
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1. For a solution to exist it is necessary to fulfil 

^32 + 1^36 — C?33 — d46 — ^44 — C^39 + (^37 + ^41 = 
^32 + ^36 — C^33 — 2(i44 + ^40 — ^39 + ^37 — dsi + ^29 + ^43 — d^e 
^32 — ^33 — 2(i44 + d^o — C^39 + ^^37 + C^SS — d^l + (^29 + '^42 — ^45 
^36 + 2^32 - dss - 2du + d^o - ^39 + ^37 - ^34 = 
d42 + 2(i29 - C^30 + ^40 " 2^44 - ^39 + (isy - ^31 = 
^38 + dio — ^39 — ^44 — ^31 + ^29 + ^^42 — ^45 = 0. 

2. rank{ME) = 12. 

Prom 2. we get the information that the representation of dQTi\TypeE as a diver- 
gence is not unique. But the important results are the equations ( [7.6[) , because we 



obtain relations among the coupling parameters if we use ( [7.4[ ): 

07 - 2as - 09 + ai2 - ^16 + bn - 618 = (7.7) 
-2a8 - 09 - aio + an + ai2 + 614 - ^15 + hi - hs + &20 = (7.8) 
-2a8 - ag - aio + an + ai2 + 614 + ^17 + ho - 621 = (7.9) 
3 1 

-2a8 - -ag - -aio + ai2 + 613 + hr - 618 = (7.10) 
1 3 

-2a8 - -ag - -aio + an + 614 + 617 + hg + ho - 621 = (7.11) 

-aio + an + 614 + &17 + ho - hi = 0. (7.12) 
These equations are direct consequences of first order gauge invar iance. 

7.2 Divergences of Type F, G, J 



In analogy to the case of Lorentz type E we obtain linear relations among the 
coupling parameters from the types F, G and J. One finds the following 9 linear 



independent relations 

— a4 — a^ — ag — aio — ^2 + ^3 — ^4 = (7.13) 

— 2a3 — 2a6 — ag — aio — an — ai2 + ^3 — 264 + 2h — bg = (7.14) 

as — ae = (7-15) 
1 1 

-a4 + as - 2a6 - -ag - -aio - 61 - 64 = (7.16) 

-2a2 - as + ae - as + ^8 + fell = (7.17) 

— 2ai + 2a2 — a4 + 2a5 — 3a6 — aj + ag — bg — 26io + ^ii — 2612 = (7.18) 

as — ag — 6g — 612 = (7-19) 
1 3 

-2«4+ 2«5- 206 + 67 -^8-^12 = (7.20) 

-2a2 + as - ae - a8 - ^8 - ^ii = 0. (7.21) 



= 
= 

(7.6) 
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Together with the six relations from type E ( |7.7| - |7.1^ ) we get 15 hnear independent 
equations which restrict the admissible theories. By construction these equations are 
necessary for a spin-2 theory to be gauge invariant. 



7.3 Nilpotency of Q 

The gauge charge operator is by definition nilpotent (Q^ = 0). As a consequence the 
application of twice the gauge variation to every expression must vanish, i.e. 

idQfT,ix) = 0. (7.22) 

If we now use the gauge invariance of Ti to first order, we get additional constraints 
for the Q- vertex Ty-^, namely 

=0. (7.23) 

This equation gives us restrictions on the parameters of Ty-^. After a lengthy calcu- 
lation one arrives at exactly 63 linear independent coefficients. The remaining ones 
can be expressed as linear combinations of them. One might think that these linear 
dependences may produce further necessary conditions beside the fifteen above. But 
it turns out that this is not the case. As was mentioned above it is shown in appendix 
B that the arbitrariness in the Lorentz types E, F, G and J can be shifted into a form 
with vanishing divergence. In view of this one can say that the Q-vertex is unique 
modulo redundant terms. In order to deal with more simple algebraic expressions 
it is convenient to work with the full set of 98 different parameters since otherwise 
we must work with long linear combinations of the independent parameters because 
the Lorentz types get mixed. We have convinced ourself that the relation ( 7.23| ) is 
always satisfied as soon as gauge invariance to first order holds. 



8 Gauge Invariant Spin-2 Theories 

The preceding section has shown what kind of restrictions we obtain if we require 



the theory to be gauge invariant. The 15 equations (7.7- 7.21 ) we have found for the 
33 parameters ai, . . . , ai2 as well as 6i, . . . ,621 play a central role. Now we can look 
at an arbitrary solution to this set of equations. The corresponding Ti is then gauge 
invariant to first order for the following reason. We have to write the gauge variation 
of this Ti as a divergence in the sense of vector analysis. Because of the generality 
of our ansatz for the Q-vertex every term in dqTi can be uniquely identified with a 
dj-monomial in With the help of the equations from appendix A we can then 
find a unique divergence for the types A, B, C, D, H and K. For the other types we 
can also find a divergence which is unique modulo redundant terms. 
Summing up we have proven the following proposition 
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Proposition 1. Let Ti and T^-^ be given as above, furthermore let f be the following 
mapping 

f : (Vectorspace of Wick-monomials) — > [Vectorspace of coefficients ai, bj) 
aih^t^K^h + ... + b2iu>'ilif^hPp^ ^ (ai, . . . , ai2, 61, . . . , 621)- 



Let V G M?'-^ be the space of solutions to ^ - \7.2J\ ). V is an 18- dimensional subspace 



of M^^, which is characterized through the following injective linear mapping L : 



0-6,0,12, &3; ^4, &5i ^61 bj, bio, bii,bi2, bis, bu, biQ, bn, bis, big, 620, &21) ' — ^ 

h - 610 - ^ (&16 - bn + bis) , ^ (&13 + bn - bis - big) , -oq + ^63 -64 + 65 
- ^ (&6 + 3 613 - 614) - bn + (3 618 - 619 + 620 - 621) , -ae + 2 (67 + 611 - 612) , 

ae, oq, 616 - 617 + bis,]^ (-613 - 617 + 618 + 619) , ai2 + 613 + 617 - 618 - 619, -ai2 

+ 613 + 617 - 618 + 619, -ai2 + 613 - 614 - 618 + 619 - 620 + 621, ai2, -64 - 2 (67 
+ 611 - 612) - 613 - 617 + 618, 63-64-2 (67 + 611 - 612 + 613 + 617 - 613) , 

63, 64, 65, 66, 67, —611, —612, 610, 611, 612, 613, 614, — 618 + 621, 616, 617, 618, 619, 620, 621] 



Then we have the two equivalent statements: 

{Al) dQTi(x)=aX/i(^) o,nd dg(a^f{)^(x)) =0 
{A2) f{Ti)eV = im{L) 

where im{L) means the image of the linear mapping L. 

This proposition determines all gauge invariant spin-2 theories up to first order of 
perturbation theory. Among them there is the trilinear couphng in the expansion of 
the Einstein-Hilbert Lagrangian. It is given by L(0, 1, —1, —1, 0, . . . , 0, 1, 0, 0) which 

11 times 

can be written explicitly 



rph,EH 



i /i^^/i,^/^,, + i h^^^hfhZ^ + h^'^h^^h^J^' 



i.l) 



The ghost coupling turns out to be the one first suggested by Kugo and Ojima |31 
namely 



T^"'^^ = Kluiufph"^ - uPyftX^ - uPpulh^"" + . (8.2) 

From the viewpoint of gauge properties of a quantized tensor field we have obtained 
a set of 18 linear independent gauge theories. We claim that the most general gauge 
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invariant theory which is given by a hnear combination of these 18 hnear independent 
ones is physically equivalent (in the sense explained below) to the trilinear coupling 
of Einstein-Hilbert ( ^.1] ) plus the ghost coupling of Kugo-Ojima ( ^.2| ) up to first 
order of perturbation theory. 

Let Pphys be the projection from the whole Fock-space J- onto the physical 
subspace Tphys, which can be expressed in terms of the kernel and the range of the 
gauge charge operator Q by 

J^phys = kerQ/ranQ (8.3) 



(see e.g. [30, |2^). Then two S'-matrices S", S' describe the same physics if all matrix 



elements between physical states agree in the adiabatic limit g ^ 1, i.e. 

\iui{<t), PphysS{g)Pphysip) = yun{(j),PphysS'{g)Pphys'ip), V(/>,-0 G JF. (8.4) 

For theories with massless fields the existence of the adiabatic limit is a problem. 
To avoid this we work with a perturbative version of (^.4|): 

PphysTnPphys " PphysTnPphys = divergences. (8.5) 

Obviously (^^) for all n implies ( |8.4| ) if the adiabatic limit exists. Specializing to 
first order n = 1 we see that two couplings Ti and T[ which differ by a divergence 
are physically equivalent to first order. Furthermore, if they differ by a coboundary, 
i.e. a term 

Tf = dgX (8.6) 

where X has ghostnumber ng{X) = —1, they are also equivalent because of the 
equation 

PphysidiQ X^ Ppfiyg PphysQX Pphys -\- PphysX Q Pphys (^•''') 

since by inspection of ( ^.31) we have 

QPphys = = PphysQ- (8-8) 

Let us return to the space of solutions V from proposition 1. Every vector in 
V corresponds through the mapping to a gauge invariant theory to first order 
of perturbation theory. As was mentioned earlier the trilinear coupling of Einstein- 
Hilbert lies in this space. Now we look at the other theories beside the Einstein- 
Hilbert coupling. For this purpose we choose a suitable basis in V. It turns out 
that a basis can be choosen so that all theories beside the classical Einstein-Hilbert 
coupling consists of divergences and coboundaries only, i.e. we have the following 
theorem: 

Theorem 1. Up to first order of perturbation theory the most general gauge invari- 
ant trilinear self-coupling of a quantized tensor field h^^[x) is physically equivalent 
(in the sense described above) to the one obtained from the expansion of the Einstein- 



Hilbert Lagrangian (given by ( \8. J\ ) without the two divergence terms, see p^/j. 
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The proof of this theorem is as follows: With the notation of proposition 1 we 
choose a basis {vi, . . . , vn, Veh) in V which displays the vector 



Veh ■■= L{0, 1, -1, -1, 0^_^, 1, 0,0) eV 

11 times 



(8.9) 



corresponding to the Einstein-Hilbert coupling with Kugo-Ojima ghost coupling ex- 
plicitly. We can choose the remaining basis vectors vi, . . . , vn in such a way that 
they have the following property: 

f~^{vi) = ^ dgX + divergences, Vi = 1, . . . ,17 

where X is of the form 

X ~ 5 I uhh or X ^ d \ uuu. 
We consider the following set of vectors {vi \ i = 1, . . . , 17} e V : 
i;i = (0,0,-l,-l,l,l,0^_^) 

27 times 

V2 = (0^, 1, -1, -1, 1, 0^_^) 

8 times 21 times 

vs = (0^_^, 1, 1, 0, 0, 1, 0^_^) 

13 times 15 times 

= (0^_^, -1, -1, 0, 1, 1, 0^_^) 

12 times 16 times 

i;5 = (0^_^,l,2,0^_^) 

16 times 15 times 

V6 = (0, 0, 1, 0^_^, 1, 0^_^) 

13 times 16 times 

V7 = (0^_^, 1, 1, 0, -1, -1, 0^_^) 

18 times 10 times 

vs = (0, 0, 0, 2, 0^_^^, -2, -2, 0, 0, 0, 0, 1, 0, 0, -1, 0^_^) 

8 times 11 times 

i;9 = (-l,0^_^,l,0^_^) 

20 times 11 times 

vw = (0^_^, 1, 0, -1, -1, 0, 1, 0^_^) 

18 times 9 times 

vn = (0,0,-3/2,0,0,0,0,0,1,1,1,1,0,-2,-1,-1,0^^^^^,1,0,0,0,0,0,1,0 

8 times 

vi2 = (0^_^, 1,1, 0,0, 0,0, 0,1) 

25 times 
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7;i3 = (0^_^,1, 0,0, 0,0, 0,-1,0) (8.24) 

25 times 

vii = (0, 0, 1/2, 0, 0, 0, 0, 0, -1, 1, 0, -1, 0^_^^, 1, 0) (8.25) 

19 times 

vi5 = (0^_^, -1,0, 1, 1, 0, -1,0, 0, 0) (8.26) 

24 times 

vi^ = (-1/2, 0, 0, 0, 0, 0, 1, 0^_^, 1,0, 0, 0, 0, 0) (8.27) 

20 times 

t;i7 = (0^_^, 1,0, -1,-1, 0,1,0). (8.28) 

26 times 

It's easy to see that these vectors together with Veh form a basis of V . What remains 
to be done is to show that they indeed have the property ( |8.10D . After a lenghty 
calculation we have found: 



f-\v^) = - uP^u^h^'' - uiu^hP^ + uPpU'ihP'' + uPpuPhPT 

= d,,{uPpUphP^ - u'^pUphPP) 
r\v^) = uPpU^K' + 2uPuPphP;^ 

= i dQ (n'^nr^u'^) + {uPpUph^"' + uPuPh'"' 
+ uPu^'hP^ - urpUphPP - u^uPphPP) 

f-\v(^) = hi^h^ph+up^u^hpr 

= -idQ{uPhPi'h) 

rHvj) = uP^uP^h + uP^uPh,^ - uPpuP^h - uP^uPh,^ 

= ^ ^^^.hh + uPh,^h) + hPJ'hh + uPuPhj 

r\vs) = 2Kj^hPPh,^ -2uiuP,hP^ - 2uP^uPhP;' + uP^uPh - uP^uP^h 
= idQ(^ ut^hh + i uPh,^h + + uPhP^h 

- uPhP^h,^^ + df, hP^hh + hP^h^^h + uPpv!'h!'P 

- 2 uPuPhP" + - uPuPh - - uPuPh + - uPuPh 

.V ' 2 P 2 P 2 



8.29) 
8.30) 
8.31) 
8.32) 

8.33) 

8.34) 
8.35) 



(8.36) 
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f-Hv9) = -h^;'h„h + uppuf^h 

i . .-„.,x 1 . .,„„,,^ (8-37) 



dQ{uf^hh)--d^{h^J'hh) 

uPuth... 

(8.38) 



d^[uPuPph - uPuPph) 
3 
2 

uluPhPp^ - uPu\ 
2 uP^hP'^h'"' + ^ uPhP^h - i uPurpU") + dc, (hP''h':^hPf^ 

+ hP^h'^'^hPl^ + ^ uPpuPhP"^ + ^ u''£t(;/i'^" - ^ ^i^n'^/if;" - ^ uZuPhP" 

- \ u^uPhP'' - ^ uP^uPhPP - ^ uPuP^hPP + ^ uPuPhP^) 



^.39) 
^.40) 

^.41) 



r\vi2) = uPpuPhi^ + uPuPphi^ + n'^^/if; 
= -idQ{uP,uPu'') 

rHvu) = \ hfh^pPh - hP^K^hP^ + - hP^h^^hP^ + <n^/if; 

= i dq (uPhP'^h'^J - ^ uPhPJh + ^ uPuPu''^ (8-42) 

- dc,[hP''h':^hP^ - hP''h'"^hPf) 
f-'M = - uPpuPhP'' + uPuPpK^ + uiuPphP'' - uPuPhP^^ 
= dp{uPulhPP - uPuPhP^) 

f-'M =-\ hP;h„h + hP^hfh^^ + uiuPphP^ 

= . (i uPph^^h'^^ - i uPphh) +dp[\ hPTh'^^h'-f' - i hPJhh 

r\vi7) = uPuPpK^ - uiuPhP;: - uPuPhP^ + uPuPhPp'' 

= ido(- uPhh + - uPh,„h - - uPh'^^h''^ - uPh'^^h'^^ + uPuPu" 

+ dp hP^hh - ^ hPJ'h'"^h"^ - uPuPhP"" + uP^v!'hPP + uPulhPP 

(8.45) 

This shows exphcitly that indeed all the basis vectors Vi are a sum of divergences 
plus coboundaries. It should be noted that there is no possibility to write Veh in the 
form ( |8.1[l| ). Then the theorem is proven because all basis vectors except Veh have a 
form which lead to trivial 5'-matrix elements between physical states. Together with 
the discussion preceeding the theorem we claim that the only physically relevant 



^.43) 
^.44) 
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theory is the couphng of Einstein-Hilbert. 
Now there arise two questions: 

Question 1: Will the statement of this theorem remains true in higher orders? 
To answer this question we have to show that in each order n we can achieve the 
form 

Tn = T^" + dq {Xn) + divergences (8.46) 

where will be constructed from , j = 1, . . . , n — 1 only. We are quite sure 
that this is indeed the case so that the divergence or coboundary contributions will 
have no physical effect. 

Question 2: What about the gauge invariance of the Einstein-Hilbert coupling in 
higher orders? In |39| Schorn obtained the result that the Einstein-Hilbert coupling 



in combination with the Kugo-Ojima-coupling for the ghosts is gauge invariant to 
second order, see sect. 5.7 and 5.8 of [37|. There it was necessary to introduce 



normalization terms which coincide with the four graviton coupling obtained from 
the expansion of the Einstein-Hilbert Lagrangian. Higher than second order have 
not been investigated up to now. 

To summarize, we have given a detailed analysis of the gauge properties of a 
quantized tensor field. Very strong restrictions on the admissible form of the inter- 
action are obtained through the requirement of perturbative gauge invariance even 
to first order of perturbation theory. Among all solutions to our set of equations only 
the Einstein-Hilbert coupling remains as a physically relevant theory. This fact is 
very remarkable since in our approach only the gauge properties of a quantum field 
describing a spin-2 particle were considered and no use was made of any geometrical 
input from classical general relativity. In view of this and with the preceding work 
about Yang-Mills theories in mind we have seen that the principle of operator gauge 
invariance is really universal. 



9 Infrared Behaviour of Massive Scalar Matter Coupled 
to Gravity 

Now we come to the second part of this work in which we investigate the long 
range behaviour of massive scalar matter coupled to the quantized gravitational 
field. The study of the infrared problem in quantum gravity goes mainly back to 
the work of Weinberg who has investigated the infrared behaviour of virtual 
and real soft graviton emission processes. The transition amplitude for a process 
in which an emission of a single graviton occurs is logarithmically divergent. In 
the sum over an infinite number of soft gravitons emited he claims that infrared 
divergencies in the transition amplitude cancel to every order of perturbation theory. 
In his treatment he considers virtual and real bremsstrahlung seperately. For the 
virtual ones he obtained in the sum over an infinte number of bremsstrahlungs 
processes a transition amplitude which is proportional to some power of the lower 
cutoff parameter A. So this rate will vanish in the limit A — > 0. For an infinite 



33 



number of real processes he obtained a transition ampHtude which is proportional 
to the same negative power in A. So the cutoff disappears if virtual and real processes 
are combined. From our point of view we think that the problem cannot be treated 
in this rather simple way for the following reasons. The structure of the problem is 
the same as in quantum electrodynamics (QED). There it is well known that the 
infrared divergencies cancel in lowest order but in higher orders one has to deal with 
sub divergences in an arbitrary Feynman diagram. The question of wheather or not 
these divergencies cancel to every order of perturbation theory is a topic which had 
long been investigated, see ||5l|, |2l|, but is still under discussion [^]. As far as we 
know the cancellation of these divergencies is not at all clear. So we think that one 
has to investigate these processes carefully order by order. 

We will reinvestigate the infrared problem in quantum gravity in the framework 
of causal perturbation theory. Especially we consider the emission of a real soft 
graviton in two particle scattering. We use the method of adiabatic switching with a 
test function from the Schwartz space which cuts off the interaction at large distances 
in spacetime. This makes all expressions well defined during the calculation and it 
turns out to be a natural infrared cutoff. We therefore avoid the introduction of a 
graviton mass. The latter, although widely used in the literature, is less satisfactory 
on physical grounds because it modifies the interaction at short distances, too. It 
is not at all clear that the different infrared regularizations give the same result for 
observable quantities. To study this question we want to calculate the differential 
cross section for bremsstrahlung from first principles. 

In order to set our notation let us start with a very short review of the quan- 
tization of the massive scalar field (p{x). The field (p{x) obeys the Klein-Gordan 
equation 

{a + m'^)ip{x) = (9.1) 
and it has the following representation in terms of annihilation- and creation-operators: 

(/9(x) = (27r)"^/^ / f' ^ {a{p) ex.p{-ipx) + d{p)'< ex.p{+ipx)) (9.2) 
J V2u;(p) 



The energy a;(p) is as usual given by ojip) = \/^~-\-rn^ . The operator a{p) an- 
nihilates a particle with momentum p whereas a(p)^ creates an antiparticle with 
momentum p. These operators obey the following commutation relations 



[a{p),a{pi]=6^'\p-p') 
[d{p),d{p')^]=6^^\p-p') 



(9.3) 



and all the other commutators vanish. The energy-momentum tensor of (p{x) is given 

by 



Tifl'ix) = : (^(p{x)'l^ip{x),^ + (p{x)l<f{x),^ - rjf,^{^{x)]^^{x),^ - m^ipix)^ ^{x))^ : 



(9.4) 
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The interaction with the gravitational field will be described through the following 
first order coupling |24| 



ri(x) = : /i"^(x)6„;3^,rr(x) 



(9.5) 



where n = \/327rG with Newton's constant G. This leads to the following explicit 
form of the first order interaction 



Ti{x) = \^\^- /i"^(a;)(p(x)L(p(x),;3 : -m^ : h{x)^{x)^ ^{x) : 



(9.6) 



This first order coupling will be used in the following subsections for the calculation 
of cross-sections. 



9.1 Bremsstrahlung 

For the calculation of the bremsstrahlung process we need to know the explicit form 
of the S-matrix up to third order in the coupling constant k. This means that we 
have to construct the time-ordered product T3(xi, X2, X3). According to the inductive 
construction of Epstein and Glaser we proceed by first calculating T2(xi,X2). We 
start with the first order term( |9^ ) and apply the inductive construction as described 
in section 2. First of all we have to build the two distributions 

A'^{xi,X2) = -Ti{xi)Ti{x2) (9.7) 

and 

R'2{xi,X2) = -Ti{x2)Ti{xi). (9.8) 

Since A'2 and R'2 are products of normally ordered field operators we apply Wick's 
theorem to obtain a normally ordered expression. We are interested in scattering 
processes so we have to concentrate on terms with zero or one contraction only. The 
contraction of the free field operators are defined by 



/1"^(X1)/1'^^(X2) 
I 1 

ip{xi)ip{x2) 

where {x) is the positive frequency part of the massive Pauli- Jordan distribu- 
tion. For the splitting operation we need the difference 

D2{XI,X2) = R'2{XI,X2) - A'2{XI,X2). (9.10) 



ib^l^l^'^ D^^\xi - X2) 
iZ)(+)(xi-X2) 



(9.9) 
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This distribution has causal support and it is given by 



+ 4 : h''^{x2)hP''{xiMx2)lip{xiU ■■ dfd;^Dm{xi - X2) 

- 2rr? : h''\x2)h{x{)^{x2),^^{x{)^ : O^'Dmixi - X2) 

h''^{x2)h{xi)ip{x2)lv{xi) : dfDm{xi - X2) 



2m' 
2m^ 



hix2)hP''ixiMx2Mxi)l : d^'Dmixi - X2) 
h{x2)h^{xi)if{x2)^ip{xi)^ : d^.'Dmixi - X2) 



+ : h{x2)h{xi)(p{x2)(pixi)^ : Dm{xi - X2) 
+ : h{x2)h{xi)ip{x2)^ ip{xi) : Dm{xi - X2) 
+ 4 : ip{x2)]pV{x2)Mxi)^^v{xi),p : b'^'^P" Do{xi - X2) 



+ 2m^ 
+ 2rr? 



- 4m 



ip{x2)]p(p{x2),p(p{x{f ip{xi) : Dci{xx - X2) 

ip{X2)^^{X2)^{xi)\^ip{xi)n : Dq{xi - X2) 

</?(x2)V(a^2)¥'(a;i)V(a^i) : Dq{xi - X2 



(9.11) 



To obtain the time-ordered product T2 we have to spht the numerical distribution 
in D2 according to it's singular order. The singular order of the Pauli-Jordan distri- 
bution is a; = —2, so the splitting is trivial and we can use the formula 



Drn{x) = DZ\x)-D-^{x). 



(9.12) 



The retarded distribution R2{x\,X2) is then given by 1)2(3^1,3^2) if we replace all 
Pauli-Jordan distributions by the retarded ones Z)''^*. Finally the time-ordered dis- 
tribution T2{x\,X2) is given by 



T2{xi,X2) = R2{Xi,X2) - R2{xi,X2). 



(9.13) 
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For later reference we here give the exphcit expression 



T2(xi,X2) = [4 : h''^{x2)hf%xi)^p{x2)M^i)]p : dl^d^^D^{xi - X2) 
+ 4 : h-''{x2)h'"'{x,Mx2)l^{xi)^ : dfdf^'Dlix^ - X2) 

- In? ■ h''(^{x2MxiMx2)M^l)^ ^a'^m(^l " ^2) 

- 2rr? : /i"^(ar2)/i(xi)(/.(x2)L<^(xi) : dfD^{xi - X2) 

- 2m2 : h{x2W{xiMx2Mxi)l : d^'D^{xi - X2) 

- 2m2 : h{x2)hP''{xiMx2)^ip{xi)^ : d^^W^ixi - X2) 
+ : h{x2)h{xi)(p{x2)^{xi)^ : D^{xi - X2) 

+ : h{x2)h{xi)(p{x2)^ ^{xi) : D!^{xi - X2) 

+ 4 : ^{x2%v{x2)Mxi)lv{xi),p ■ r'^P'^D^ (xi - X2) 



(9.14) 



+ 2rr? : Lp{x2)]p^p{x2),p^p{xi)^(p{x{) : D^{xi - X2) 
+ 2iT? : ip{x2)^(p{x2)ip{xi)l^ip{xi)^ : D^{xi - X2) 

- 4m^ : ip{x2)^ (p{x2)ip{xi)^ ip{xi) : Dq{xi - X2) 

- 4 : h'''^{x^)^{xi)l^{x^),phP%X2Mx2%v{x2)^ : 
+ 2m^ : h°'l^{xi)ip{xi)]^^p{xi),i3h{x2)(f{x2)^ip{x2) : 
+ 2rr? : h{xi)ip{xi)^ ip{x-i)h'"' {x2)ip{x2)^,p(p{x2)^ : 

- : h{xi)ip{xi)^(p{xi)h{x2)(p{x2)^fix2) : 

To obtain the time-ordered product to third order we proceed in much the same 
way as in the calculation of T2. We calculate the distributions R'^{xi,X2,X3) and 
A'^{xi,X2,Xs) which are given by 

R3{xi,X2,X3) = T2{xi,X3)Ti{x2) + T2 (x2, X3)Ti (xi) + Ti (a;3)T2 (a;i , X2) (9.15) 

and 

A'^{XI,X2,X3) = fi{x2)T2{xi,X3) + fi{xi)T2{x2,X3)+T2{xi,X2)Ti{x3) (9.16) 

where Ti and T2 are defined by 

Ti{xi) := -Ti{xi), i = 1,2,3 (9.17) 

and 

T2{xi,Xj) := -T2{xi,Xj) + Ti{xi)Ti{xj) + Ti{xj)Ti{xi), i,j = 1,2,3. (9.18) 

It should be noticed that we only collect those terms which have two contractions 
involving all the spacctimc points xi,X2 and X3, i.e. there is no product of contraction 
functions with the same argument. After we have build the causal difference D3 = 
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i?3 — ^3 we split the numerical part into Rs with support in the backward light cone 
and ^3 with support in the forward light cone. The time-ordered distribution is 
then given by 



T3{xi,X2,Xs) = R3{xi,X2,X3) - R'^{Xi, X2, X3) . 



(9.19) 



Because the calculation is staightforward and not very enlightning we only give the 
result 



IK, 



T3{xi,X2,X3) = ^ ^ 



§[: /i'"'(2:7r(l))¥'(a;7r(l))t¥'(a;7r(2))iy¥'(a;7r(2)).£¥'(a;7r(3))« : 



+ : ^'"^(a^7r(i))v'(2;7r(i)).aV5(a^7r(2))V(^'r(2))'£V'(a;7r(3))L : 
X 6"^^^I)o''^(x,(2) -^.(3))5;^'''5p-'^'^^(^4i) -a:.(3)) 
+ 4m^ : /t'"'(a;^(i))(/c(x^(i))J,(/?(a;^(2))V(a;7r(2))</5(a;7r(3))'a : 
+ : ^'"'(a;7r(i))¥'(x7r(i)),p¥'(a:;^(2))V(a;7r(2))¥'(a^7r(3))L : 

X Dq{x,^(^2) - a;,r(3))9a"*^^<9a"*''i^m(a;7r(l) " 2:^(3)) 

+ W : h''" {x^(i))ip{x^(^i))l(p{x^(^2))l^{Xn{2))'e^{x^{3)) ■ 

+ ■ h'"'ix,^{l))(pix^{l)),p(pix^{2))i'PixTr{2))'e'P{x^3))'' ■ 

X -D^(a;7r(3) - a;,r(2))5a''*'^£>m(a^7r(i) - a;,r(3)) 

+ W : /i(a;^(i))<^(x7r(i))V(a;7r(2))V(^'r(2)).£'^(a;7r(3))'/3 : 

+ : Hx^^i))ip{x,„(^^))tf{x^(^2))ttfiXni2))^f{Xn(3))}p ■■ 

X 6"^^''Z)^(X^(3) - X^(2))<9q''*'^L'^(x^(i) - .T^(3)) 

+ 2m^ : /i(a;,r(l))¥'(2:7r(l))V(2;7r(2))V(^T(2))r7</^(a;7r(3)) : 

+ : h{x^(l))f{XT,(l))<pix^^2))ly<Pix^{2))-r'Pi^7:{3))^ ■ 

X Dq{x^^3) ~ .X^(2))-Dm(x^(i) - ,T^(3)) 

+ 2m^ : /i(a;^(i))v?(a;,r(i))V(a;7r(2))V(2;7r(2))9'(a;7r(3))'/3 : 

+ : h{x7r{l))'P{Xn{l))'Pix^{2))''^{Xn{2))vixT,(3))]p : 
X ^o'(a^7r(3) - 2:7r(2))^^"*^^^m(a^7r(l) " 2:7r(3)) 

+ 8m^ : ^'"'(a:^(i))95(-^^7r(i))fp95(-^7r(2))V(a;7r(2))</'(a;^(3)) : 
+ : hf"'{x,„(^i))<p{x^i)),p<f{x^^2))^ip{x,„^2))'P{x^{3))'' ■ 
X ^o'(3^7r(3) - 2:7r(2))9a"'''-D^(a;^(l) - a;^(3)) 
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+ 4m^ : /i(2;^(i))93(a;^(i))V(2;7r(2))V(2;vr{2))¥'(a;,r(3)) : 

+ : KXn(l))'P{x^(i))ip{x^(2))^Lp{x^(2))v{Xn{?.))^ ■ 



(9.20) 



Due to the sum over all permutations of the indices this expression for T3 is 
obviously symmetric in it's arguments as it is required by the definition of the time- 



ordered product 1 29]. The third order 5-matrix for the bremsstrahlungs processes is 
then given by 



^3(5) 



d'^xi . . . d'^xsT3{xi,X2,X3)g{xi)g{x2)g{x3), g e S{R'^). 



(9.21) 



Now we consider the scattering process of two massive scalar particles in which 
a bremsstrahlungs graviton is emited in the final state. That is, we want to calculate 
the expectation value of 5*3 ((7) between the following initial and final states 



d^pid^qiMPi)MQ'iHpi)^a{qi)^n (9.22) 
d3fcd3p2d^g2^(p2, g2)0''"(^)a^"(^)^a(p2)^a(g2)^^^. (9.23) 



These are wave packets in momentum space where ipi, i = 1,2 are one particle 



wave functions from L' 



^' is a two particle wave function from 



and (f)^'^ is a tensor- valued square integrable wave function on M . The vector Q. is 
the vacuum vector in the asymptotic Fock-Hilbert space. We want to calculate the 
following S'-matrix element: 



(9.24) 



To show the details of the calculation we restrict ourself to the following term from 
the above computed T3 ( 9.20| ) 



rj^^(xi,X2,X3) := itC" : /i'"^(xi)(/?(x3),/3(/?(a;i)J,(/7(x2)J^V'(a;2).£ : 
h^^^'D^[x2 - x^)d^^^dl'D^[xi - X3) 



(9.25) 



where the (1) refers to this particular term of T3. The S-matrix element correspond- 
ing to this contribution is then given by 



S 



(1) 



d'pid'q{d-'kd'p2d'q2d''xi . . . d^X3MPl)MQl)'^iP2, {k) 

X g{xi)...g{x3){a^'''{k)^a{p2)^a{0n,T^^\xi,X2,X3)aipi)^aiqi)^n). 

(9.26) 
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Then we have to compute the vacuum expectation value 



{n,a{q'2)a{p-'2)af'''{k) : h^'^ {xi)ip{xs)M^i)^M^2)ly^{x2h ■ a{p-{)^ a{q-{)^ Q) . (9.27) 

This can easily be evaluated by inserting the Fourier representations of the free field 
operators ( |4.4| ) and ( ]9.2| ) and the use of the commutation relations in momentum 
space (^^) , • The result is given by 

exp [+i{kxi -P1X3 + q2Xi + P2X2 - qiX2)\ pf P2 I1 
+ exp [+i{kxi - qiX3 + q2Xi + P2X2 - P1X2)] q{ ^2 P2 Pi (9.28) 
+ exp [+i{kxi - P1X3 + P2X1 + ^22^2 - qiX2)] Pi P2 qI 1i 
+ exp [+i{kxi - qixs + ^22:1 + ^22:^2 - P1X2)] gf P2 qJ Pi ■ 



The S'-matrix element then reads 



Sf^ = i(27r)-l5/2^3 J d^picfq-{d'k(fp2d^q2d^xi . . . d^XsMpi)Mqi)^{p2, ^2) 

r''{k)g{xi) . ..g{x3)b''^""D^{x2 - X3)d^^^ d^' D^{xi - x^^P"" 

{32uj{k)uj{pi)uj{p2)uj{qi)u;{q2)y^^^ 

exp[+«(A;j;i - ^1X3 + 922;! +P22;2 - gi2;2)] Piq2P2ll 

+ exp[+i(/ca;i - 51X3 + Q'2a;i +^22:2 -Pia;2)] q{q2P2Pi 

+ exp[+i(A;xi -^1X3 +p2a:;i + ^22^2 - gi2;2)] Pi P2 ^2 9i 

+ exp[+i(A;xi - ^1X3 +P22;i + q2X2 -Pi2;2)] P2 ^2 f*! • 

(9.29) 



We introduce the abbreviations Ii, ■ ■ ■ , /4 for the four parts of S^j^ . Then we consider 
the spatial integrations in Ii. They can be carried out by inserting the Fourier 
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transforms of the test functions as well as the Feynman propagators: 

h= J d'^xi . . . d'^X3g{xi) . . . g{x^)D^{x2 - X3)d^' d^^ dI^{xi - X3) 

exp[+i{kxi -P1X3 + 922:1 +P2X2 - qiX2)] 

= (27r)-^'^ j d^xi . . . d^X3d\ . . . d%g{h) . . . g{h) exp [-i{hxi + I2X2 + hx^)] 

^^(/4)exp[-i^4(x2 -X3)]5^3a^i£)^(l5)exp[-i^5(xi -X3)] 
exp[+i(A;xi — P1X3 + q2Xi +P2X2 — qiX2)\ 

= (27r)-i° j d'^xi . . . d^xsdHi . . . d%g{h) . . . gik) exp[-i{hxi + I2X2 + ^s^^a)] 

exp[-i/4(x2 - X3)]l'il'^D^{l5)exp[-ik{xi - X3)] 
exp[+i(fca;i — P1X3 + q2Xi +P2X2 — qiX2)\ 

= (27r)2 J d\ . . . d%gih) . . . g{l3)b^{k)mDl{h)5^^^ (/, + _ - q^) 

8^^^ {I2 + k-p2 + qi)S^^^ ih -U-k+Pl). 

(9.30) 

We can do the integrations w.r.t. Z4 and Z5 immediately and get 

h = (27r)2 J d% . . . d%g{h) . . . g{l3)DE {h + h + Pi - k - 52) 

[I? -p'i + k- + q^] [-11 -pl + k'' + q^]D^{-h + k + q2) ^^-^^^ 

(5^^^ (Zi + h + h + Pi + qi - q2 - P2 - k). 



Then the first term of S^^^ becomes 



45'''^ = (27r)-i^/2.^3y Sp^d^qid^kd^p2d^q2^^{p^)MQi)^W^ 

j d% . . . d%g{h) . . . gihW^^'b^ih +h+Pi-k-q2) 

-p't + k^ + [-11 -pl + k'^ + q!^]b^A~h +k + q2) 
6'^'^''>fg2'P29f(32a;(^)a;(#i)a;(#2)a;(gi)a;(r2))"'/' 

5^^'> {h + l2 + l3 + pi + qi-q2-P2 - k) 



/ 



d% . . . dH3g{li) . . . g{l3W^^ %, ^ \ \ ^LS "'"' 

[h+k+Pi-k- q2Y + ?0 



,„2 Xh+^k + ^f~ ^^ p^,^p]gf (32u;(fc)u;(p0c^fe)c.(gi)u;(g--2))-^/^ 
(5^^^ {h + l2 + h+Pi + qi- q2-P2-k) 

(9.32) 
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where we have inserted the propagators exphcitely. In this expression we can carry 
out the adiabatic hmit in the variable I2 since it doesn't appear in the argument 
of the propagators. We do this in the fohowing way. First of all we choose a fixed 
test function go G 5(M^) with the property go{0) = 1. Then a new test function is 
defined by 

gs{x) ■■= go{ex). (9.33) 

The adiabatic limit — > 1 is then equivalent to the limit e ^ 0. The Fourier 
transform of g^ is given by 

9eip) = ^9o{^)- (9.34) 

Then we have 

hmlffof-) = (2vr)25W(p). (9.35) 



6^0 e^-^-Ve 

Using this in ( |9.32| ) the integration w.r.t. I2 gives just the factor (27r)^. We are 
interested in the S'-matrix element in the limit where only gravitons at low momenta 
(soft gravitons) are emited. So we omit all small quantities in the numerators, and we 
take only the leading terms in the denominators, neglecting terms o(e). This is known 
as the eikonal approximation in the literature, see e.g. j^. With this approximation 
we obtain for the argument of the massless propagator 

ipi-q2-k + h+ hf = {pi - q2f + i-k + h + hf - 2{pi - q2){k - h - I3) 

= (Pi - + o{e). 

(9.36) 

where £ is the scaling parameter from the adiabatic limit. The argument of the 
massive propagator becomes 



- (-/i + k + q2f = m^ -ql- {-h + kf + 2q2{h - k) 
= 2q2{h - k) + o{e^) 



(9.37) 



since the momentum q2 is on the mass shell. Then the S'-matrix element becomes 

^l"''^ = - (2vr)-15/2zK3 / dVl(i'glc^'fed¥2d3g2^l(m)^2(gi)^(p2,g2)'^'''^(fe) 



d\d%m)m)b'''''^^^2p7r^^^''''" 

(Pi - q2Y 2q2[li - k) 
Pi 92^2^1 (32a;(A;)a;(pi)u;(^2)^(gi)w(g2))""^^^ 
S^'^Kh + h+Pi + qi- q2-P2- k). 

(9.38) 
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Now the adiabatic limit in Z3 can be done as above and we arrive at the result for 
the ^-matrix element 



/ ^4/ \i^aO-ye (^2 ~ Pi) (^2 ~ Pi) ^nvpa 

Piq2Phi{^'^^(k)uj{pi)uj{p2)uj{qi)(j{q2))~^^^ 
(5^^) {h + pi + qi - q2 - P2 - k). 

(9.39) 

There remains to calculate the other three parts S^^^^^\ . . . , S^j-'^'^\ The calculations 
are similiar to the one above so we only present the results 

41'''^ =- (27r)-lV2^^3y dVlC^^9lti'^dV2C^^92V'l(Pl)^2(gl)*(p2,92)<^''"(^) 
/ n(l \h'^l^l^ 

(g? - <ll) iQ2 - Ql) ,nuna 

QU2P2Pi{^^'^ik)u{pi)u{p2)u{qi)u{q2))~^^'^ 



S^'^^ {h + qi+Pi-q2-P2-k) 

- (27r)-iV2^^3 J d^f^d^q-id^kd^f^2d%MPi)MQi)^{P2, qi)^''^) 



I 



h(U \h'^l^^^ ^"^ '-^^ h'^'^P'^ 



]3iP2929i(32w(^)w(Pi)w(P2)w(gi)'^(g2)) 
S^^^ {h +qi+Pi- q2-P2- k) 

^j;-'*) = - (27r)-l^/2 .^3 f d-Vld-\l(i'^f^¥2d'g2V'l(#l)V'2(gl)^(#2,g2)<A'^"(fe) 



/ 



{qi - P2)^ 2p2{li - ky 

Q1P2Q2P1 (32u;(fe)c<;(#i)c<;(p2)w(gi)a;(g2))~^^^ 
(J^^^ (^1 + + - ?2 - P2 - A;) . 

(9.40) 

Now we are ready to consider the differential cross section. 

We have to calculate the absolute square of the S-matrix clement We do 

this again separately for the four parts Sj\^'^^ \ . . . , S^j^'^'^K For the first part we start 
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with (9.39) and write 



[Q2 -PiY 2^2(^1 - k) 



P'\ (I2 P>'i (I'l ['i2u){k)uj{pi)uj{p2)uj{qi)uj{q2)y^''^ 
{32u;{k')oj{prMp-2'MqrMq-2')y'^^ 

ih +Pi+qi-q2-P2- k)d^^^ {I2 + p[ +q[-q'2-P2- k') 

(9.41) 



where a * indicates complex conjugation. We use the completeness relations of the 
final states, i.e. 



^{P2,q2)^{p2\ q2T = S^^^ (p^2 - p-'2)S^^^ {q2 - q2) (9.42) 

/ 

^(^^'^(A?)(^^''^'(fc')* = <''<'^'<5(^)(fc - P) (9.43) 
/ 

where the sum is taken over all final states. Then we have 

^l^l'^'f = (2^)-"^^ I d''pid''qid^kd^p2d%d^p-i'd^qrMPl)i'2{q~i) 
f 

Up-iTMqir I d\d%g{hrg{l2W^''^^^^^^^^^ 



(q^ - p^)2 2q2{h - k) 



P P 1 eva'a'-y'e'iq2' - P'y ) iq2 - P'l ) I P' p' 7' I e' 
^'^^^^^^^ {q2-p\Y 2q2{l2-kf ' ^^^^^^ 

■^^^^ {h +Pi+qi- q2-P2- k)S(^^ {h + p'l + q[ ~ q2 - P2 - k) 

\J {32uj{k)uj{pi)uj{p2)uj{qi)Lu{q2)) {32uj{k)uj{pi')uj{p2)Lu{q[')uj{q2)) 

(9.44) 



Now we assume that the wavefunctions tpi and tp2 are sharply concentrated around 
the initial momenta pi and qi, so that we can replace the momenta in the propagators 
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by these initial momenta. Then the expression can be simphfied to 



f J J 

iQ^-p?) {^-pD Bo,..a'SWe' {Q^' ' pf) ' p{ ) 

(92 - Pr? 2q2{h - kf' (52 - Pi^ 2^2(^2 - k) 

Pi 92 P2 Qi [^2i^{k)u;{pl)uj{f2)uj{qi)uj{q2)] 

I d^pid^qid^p-{'d^q-{'i>i {pi)Mqi)i^l (pi')*V^2 {qi'T 

{h +Pi+qi-q2-P2- k)5^^^ {h + p'l + q[ -q2-P2-k). 

(9.45) 



We observe that the last integral herein depends on the initial state only. We denote 
it by F{P, hjh), where we have introduced P = P2 + q2 + k. We have 



F{PMM) = j d^pid%d^prd^qrMPi)MQi)MP~'iTMQiT 

S^'^'^ {h +Pi + qi-q2-P2- k)S^^'^ ih + p[ + q[ - q2 - P2 - k) 

= {2tt)-^ j d%dS2d^pi...d''qrMPi)MQVMP~'i'rMQ'iT 
cxp +Pi+qi- P)yi] cxp [+iil2 + p'l + q'l - P)y2] 

= J d'^yid'^y2i}iiyi)i)2iyi)i'i{y2)*i^2{y2)* ew[-ihyi]exp[+il2y2] 
exp[-z(yi - y2)P]. 

(9.46) 

The function F is normalized according to 

F{PMM)d^P = C^T^? j A|V'i(y)lV2(2/)l'exp[-i(/i -^2)y] (9.47) 



and it is concentrated around -P = P2 + 92 + ^ ~ Pi + 9i- In the limit of infintely 
sharp wave packets we may represent it by 

F(P,Zi,Z2) = (27r)2<5(4)(P-p,-g,) j A|Vi(y)|'|V'2(y)|'exp[-z(Zi-Z2)y]. (9.48) 
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This will be inserted in ( |9.45| ) and we obtain 



f J J 

{(j - P?) (12 - Pi) „/^„P„7„er,«'/3'V£' (12 -Pf') {4 -Pi) 

{q2-p^)^2q2{h-kf'^^^^^' {q2-p^)' 2q2{l2-k) 

Pi 12 pI 4 [32cj(fc)a;(pi)u;(pi)a;(gl)a;(g^2)] 

5(')(P2 + (?2 + fe-pi-gO \ AlV'i(y)l'IV'2(y)|'exp[-i(/i-/2)y]. 

(9.49) 



We can set l\ = I2 in the last integral because we neglect all contributions o(e) in 
the numerator, see the remarks after ( |9.32| ). We denote the absolute square of Sfi 
by Vfi •= I'S'/iP- Then the definition of the cross section is as follows: We consider 
a beam of incoming particles of radius K incident on the target which is assumed 
to be at rest. The transition amplitude p/j must then be averaged in space over the 
cylinder of the beam of incoming particles. Then the cross section is given by 

a = Jim 7ri?2 (9.50) 



where we sum over a complete set of final states. Then we restrict ourselves again 
to ( |9.25 ) of T3 and furthermore to the first integral 1\ ( |9.30D of If we denote 



the integration variables p2 and g2 by and gj, respectively, then we obtain for the 
cross section 

a = (27r)-9^6[4^(^-,)^(g-,)]-i^( gV(9l) f d%d%9{h)g{hT 

y {{PiQi)'^ - m^) 
f d^k d\-f d^p-f . {qj-pt) {q'f-pD s 

J 2u{k) 2u:{q-)) 2a;(p» {qf-p^)^ 2qf{h - kf^ ^'^^ (9.51) 
{qf -P?') {q'^f -Pi') 13' p' ^' e' ,(3) ^ ^ ^ r - -^ 

S(Lo{p~f) + Lj{qf) + uj{k) - uj{pl) - uj{ql)). 



To carry out the remainig integrations we have to know the integrand explicitely. 
So at this point of the calculation the tensor structure of the integrand becomes 
important and we have to write down the various contractions between the four- 
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vectors pi,qi,pf,qf. With the tensor b^^'^'"^ (4:-3) we get for the cross section 



9^6 



- m 







d%d%g{h)g{h 



(fk d^q) d^pf 
2uj{k) 2^(9/) 2u;(p7) 



{uj{qf)uj{p~f) - qfpf) {piqi) + (u;(g7)w(gl) - g/gj) (u;(pi)a;(p7) - piPf) 

- {^{(ff)^{Pi) - QfPi){^{Pf)^{Qi) -PfOi) - {^{Pi)^{Pf) -PiPf){pm) 

- {pm){^{Pi)^{pf) -PiPf) +Pi{^{pf)^{qi) -PfOi) i^ioff - q/) 

- {^{<lf)^{Pf) - QfPf) (PiQi) - {^{<lf)^{qi) - QfOi) {i^{Pi)^{P~f) - PiPf) 
+ [^{<lf)^{Pi) - QfPi){^{Pf)^{Qi) -PfQi) + {y^{Pi)^{pf) -PiPf){PiQi) 

+ {piQi){^{pi)^{p~f) - PiPf) - p}{^{pf)^{qi) - Pf<ii) i^{Qf)^{pi) - (ffPi) 



S'-^Hpf + qf + k- pi- qi)6{uj{p~f) +oj{qf) + uj{k) - uj{pi) - uj{qi)) 
{qf - PiY2qf{h - k){qf - PiY2qf{l2 - k) 



(9.52) 



It is convenient to choose the center-of-mass system, defined hj pi + qi = We can 
easily do the integration w.r.t. pf. 



a = 



9^6 



{{P^ 



m 







d%d%g{h)g{l2y 



d^k d^q) 1 
2uj{k) 2u;(g7) 2cj(g7 + k) 



{uj{qf)uj{qf + k)+ qjiq) + k)) {piqi) + {uj{qf)uj{ql) - q~fql) {uj{pl)uj{qf + A;) 

+ Pi{q~f + k)) - {uj{qf)uj{pi) - qfPi) {uj{qf + k)uj{ql) + {q~f + k)ql) 

- {Lj{pl)Lj{q~f + k) +Pi{q~f + k)){piqi) - {piqi){uj{pi)uj{q~f + k) + Pi{q~f + A;)) 
+ pf{uj{qf + k)uj{qi) + {0 + k)ql) {u}{qff - q/) - {u{qf)u{qf + k) 

+ (ffiqf + k)) {piqi) - {u{qf)u){qi) - (ffql) {uj{pl)uj{qf + k) + pi{qf + A;)) 
+ {^{<if)^{Pi) - QfPi) i^iOf + k)uj{qi) + (q) + k)ql) + {uj{pl)uj{qf + k) 
+ Pi{<lf + k)){piqi) + {piqi){uj{pi)uj{qf + A;) +Pi{qf + A:)) 

- Pii'^i'lf + ^)^(^i) + iOf + k)qi) {Lj{qf)u{pi) - (ffpl) 

5{ijj{qf) + oj{qf + k) + uj{k) - 2uj{pi)) 
{qf - PiY2qf{h - k){qf - Pif2qf{h - k) ' 

(9.53) 
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Now we are left with the integration w.r.t. qj. In order to be able to do this we have 
to rewrite the argument of the delta distribution showing the dependence on \qf\ 
explicitely. We want to use the formula 



n 



(9.54) 



where the Xi are simple zeros of the function /. For that purpose we define 



/(qf) = V g/ + ^2 + -^(g) + A;)2 + m2 + uj{k) - 2w(pl) 

= \/ qf2 + 771.2 _j_ qf 2 -|- k2 + 2qf k COS 6* + m2 + k — 2uj{pi) 



(9.55) 



where 9 = <{qf,k) and we have used the notation k = |A;|, qf = \qf\. In order to 
simplify the notation in the following formulas we write oj = 2u{'Pi). The zeros of 
the function / are given by 



qf«(k,0) = ^(k,0)+i?(k,0) 
qf(2)(k,0) = ^(k,^)-S(k,e) 



(9.56) 



where we have introduced 

Lo{2k — w)kcos 9 



A{k,9) 



B{k,9) 



2(cj-k)2-2k2 cos2 



(a; - k)^/uj^ - Aij^k + Sujkm'^ - 2k'^m'^ + Auj'^{\<? - m?) + 2k2m2 cos 26* 



2(a;- k)2 -2k2 cos2 



(9.57) 



Clearly these solutions are themselves functions of the momentum k and the angle 



9. If we consider the limit of these functions ( 9.56 ) as k goes to zero we obtain 

limqf(^) = ±Pi, j = 1,2 (9.58) 

where pi = \pi\. The derivative of / with respect to qf is given by 



qf 



+ 



qf + k cos 9 



t^qf ^Jm? + qf2 Y^k2 + m? + qf2 + 2kqf cos i 



(9.59) 



By inserting the two zeros of / into ( |9.59| ) we observe that, by taking the limit k — > 0, 
we get a finite result, namely 



rf/(qf) 



lim 
k^o aqf 



qf=qf 



(i) 



2pi 



(9.60) 
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Now we can rewrite the cross section as 



1 

"=4 



d\d%g{h)g{l2) 



{{piQiY - m'') J J 2uj{k) M<ff) 2u{q-f + k) 

{ijj{qf)uj{qf + k)+ qf{qf + k)) (piQi) + {(jj{qf)(jj{qi) - (ffql) {uj{pl)oj{qf + k) 

+ Pi{qf + k)) - {u{qf)u{pi) - qfpi) {uj{qf + k)u{ql) + {q} + k)qi) 

- {uj{pl)u{qf + k) +Pi{qf + k)){piqi) - {piqi){uj{pl)uj{qf + A;) +Pi{qf + k)) 
+ {^^{0 + k)u{ql) + {qf + k)qi) {u{qff - qf'^) - {u{qf)u{qf + k) 

+ QfiOf + k)){piqi) - {uj{qf)uj{q:i) ~ (ffqi) {uj{pi)uj{qf + k) +Pi{qf + k)) 
+ {^{'lf)^{Pi) - QfPi) i^iOf + k)uj{qi) + (^7 + k)ql) + {uj{pi)uj{qf + k) 
+ PiiOf + k)) (piqi) + (piqi) {uj{pi)uj{qf + k)+ Pi{qf + k)) 

- Pi{^i0 + k)uj{qi) + {qf + k)ql) {u{qf)uj{-pi) - qfPi) 

1 A 1 



{q}-Pr?2qf{h - h){qj-p,f2qf{h - k) j^^ |/'(qf)| 



^ I f'( 



qf=qf 



0)1 



.<5(qf-qf(^')). 



(9.61) 



Before we do the integration w.r.t. qf let us abbreviate by Fi(qf, k, kcos0) the ex- 
pression which comes from the tensor part in the numerator, i.e. 



Fi(qf, k, kcos 6) = 



{u{qf)u{qf + k)+ q-f{q-f + k)) {piqi) + {u{q-f)u{ql) - q-fql) 

{uj{p,i)uj{qf + k) +Pi{qf + k)) - {uj{qf)uj{p,i) - (ffpl) {u{qf + k) 
^ioi) + iOf + k)q^i} - {^{fd^iOf + k) +Pi{qf + k)){piqi) - {piqi) 
{u{pi)u;{qf + k) + Pi{qf + k)) +Pi^{(^{qf + k)u{ql) + {q} + k)q^) 
{uj{qff - q/) - ^{uj{qf)uj{qf + k) + qf{qf + k)) {piqi) - {u{qf) 

^{Qi) - Qm) {i^{Pi)i^{Qf + k)+ pi{qf + k)) + {uj{qf)uj{pl) - qfPi) 
{uj{qf + k)uj{ql) + {q) + k)ql) + {u;{pi)uj{qf + k) +Pi{qf + k)) 
{PiQi) + (Pm) {ojiPd^iQf + Piiq} + k)) - pf {(jj{qf + k) 



^iOi) + iOf + k)qi) {^{qf)^{m) - QfPi) 



(9.62) 



49 



Then the expression for the cross section can be written as 



J d%d%g{h)g{l2T J 



J {{piQi? - m'') J J Mk) MOf) 2uj{q-f + k) 



Fi(qf,k,kcos 9) 1 _ (j). 

(9.63) 



Prom this expression we can identify the differential cross section We do the 
qf-integration and obtain 



(qf(-^'))^ Fi(qf(^'),k,kcosg) 

2^((qf(^'))^ + 2y^((qf0))2 + k2 + kqfO) cos^ + m^) 

/ ^ y (9.64) 

\ m2 - 2V(qf(^))2 + m2a;(^i) + 2qf 0) pi cos a+pf) 

1 



2 [^(q7^))2T^/? - qf(j)li cos (3 - ^(qf (j))^ + m^k + qf(j)kcos 6] 

1 



2[^(qf0))2 + mH^ - qf0)l2 cos 7 - ^(qf^^'^)^ + + qf^^^kcos^] 



where we have introduced 



a = <{Qf,Pi) 

''=^®-'-* (9.65) 

7 = <iQf,h) 

li = l^d, i = 1,2. 



In the next step we omit all small quantities in the numerator, and we set k = in the 
denominators if there occur no singularities. Then the expression can be simplified 
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to 



da _1 {2tt)-^k^ Lo{pi) 



d%d%g{h)g{l2r 



d^k 



2 2 / 1 \ • 

1 

2[u;(pl)/0 - qfO)|k=oli cos/3 - a;(pl)k + qfO)|k=ol<cos0] 



2[uj{pl)l'^ - qfO)|k=ol2COS7 -(j(p;)k + qfO)|k=okcos(9] ' 



(9.66) 



Now one observes that the function Fi is independent of cos 9 and depends only on 
the initial momenta pi and qi so we can take it outside the integral. Then we get 



da _ 1 {2tt)-^k^ 
dn ~ 4 



d^k 



1 



2u;{k) 2[cu(pl)(/? - k) - qfO)|k=o(licos/?- kcos^)] 



(9.67) 



2[uj{pl){q - k) - qfW|k=o(l2COS7- kcos0)] 
where the function F2 is given by 

i^2(qf(^\=o) = ^-^i^i(qf(^'^|k=o,0,0) 



m? - 2uj{piY + 2qfO)|k=oPiCOSQ + 

(9.68) 



If we use ( 9.58 ) we can write the differential cross section in the form 



da 1 {2t:)-^k^ 



dn 16 



d%d%gih)g{l2 



F2(+Pi) 



d^k 
2uj{k) 



1 



\ijj{pi){l\ — k) — pi(li cos P — kcos 9)\ \u}{pi){l2 — k) — Pi(l2 C0S7 — kcos 6*)] 
d^k 



+ i"2(-Pi) j 



2uj{k) 



[uipDil^ - k) + pi(li cos/3 - kcos 6*)] [uj{pi){l^ - k) + Pi(l2COS7 - kcos 6*)] 

(9.69) 
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In order to do the integration w.r.t. k we use an identity due to Feynman 



1 



dx 



a,b £ 



ah Jq [ax + 6(1 — x)]2 ' ' 
We consider the first term in ( |9.69| ). Then we define the functions a and h by 

a = \uj{'pi){l\ — k) — pi(li cos /? — kcos6')] 
b = \uj{'pi){l2 — k) — pi(l2 cos 7 — kcos6')] . 

The integral becomes 

^ — r = t: /" kdkdr^k / — o 

2uj{k)ah 2 J Jq ^ax + b{l-xy 



(9.70) 



(9.71) 



— / kdk / dx 

2 







dcos 9 



(9.72) 



-1 [ax + 5(1 



X) 



i2 • 



We insert a and h and obtain 
vr / kdk j dx j 



dcos 6 



1 2 • 



[-iAj{'Pi)k + pikcos6' + {u){'Pi)l\ - Pill cos(3)x + (a;(pi)/^ - pib cos7)(l - x)]' 

(9.73) 



We set cos 9 = z and we observe that the integral w.r.t. z is of the form 

dz 2 



.1 [CZ + (i)2 d2 _ g2 

see 1^], where we've set 
c = pik 

Ci= -L0{'Pi)k + {uj{pi)l\ - Pill cos (w(pi)^2 - Pil2COS7)(l 

Then we obtain for ( plTSj) 

27r / kdk / dx ^ ^ ; : ; 



(9.74) 



(9.75) 



(9.76) 



where A{li,l2.,x) and B{li,l2,x) are given by 

A{li,l2,x) = - 2L0{pi) {oj{pi)li - Pill COsP)x + {oj{pi)l2 - Pil2COS7)(l - x) 

r T 2 

B{li,l2,x) = {uj{pi)li - Pill COS/3)2; + {uj{pi)l2 - Pil2COS7)(l - x) 



(9.77) 
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Then this integral w.r.t. k is of the form 



/ 



ydy 



1 



ciy2 + C2y + C3 2ci 



ln(ciy^ + C2y + ^3,) - ^ J 



dy 



2ci J ciy2 + C2y + C3 



(9.78) 



where the constants Cj, i = 1, 2, 3 are given by 

2 2 
Cl = = m 

C2 = ^(^1,^2, a;) 
C3 = B{li,l2,x) 



(9.79) 



In our case the discriminant A := 4ciC3 — C2 is given by 



A = + 4p: 



{u){pi)ll - Pill cos {uj{pi)ll - Pil2COS7)(l -x) 



'iujipif {uj{pl)l^ - PillCOS^)x + {uj{pl)l^ - Pil2COS7)(l 



4pi 



{uj{pi)li - Pill COS /3)x + (w(pi)/2 - Pil2COS7)(l 



X 



so we see that A < 0. Then the integral ( |9.78| ) is given by 



(9.80) 



ydy 



ciy^ + C2y + C3 2ci 



ln(ciy^ + C2y + C3) 



C2 



In 



2ciy + C2 - \/-A 
2ciy + C2 + 

(9"81) 



Here we consider gravitons up to an energy ujq which is much smaller than the energy 
of the incident particles. Since every real detector has a finite energy resolution below 
which he cannot detect any particles we have to integrate over all these contributions 
up to the value ujq. Then the final result of the k- integration is given by the following 
expression 



Jo 



kdk 



+ Ak + B 



1 

2^ 



ln(ciWo + ^"2^0 + C3) 



C2 2cia;o + C2 - 

" I 2CIUJ0 + C2 + 



ln(c3) + 



C2 



In 



C2 



C2 



(9.82) 



The second term in ( |9.69| ) can be treated in the same way. The only difference is 
the sign of pi in the definition of the functions a and b, se e (|9-71 ), which has no 
consequences for the calculation of the integrals ( 9.74 ) and (|9.78| ). 
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9.2 Adiabatic Limit 



In the preceeding subsection we've found the differential cross section for brems- 
strahlung in a scattering process of two massive scalar particles. Now we want to 
discuss the adiabatic limit. For the differential cross section we've found 



da 

dn 



1 



(2^)-9k6 



16 



{{piQiy 



d\d%g{h)9{hrF2{+9\) 



27r / dx- — 
Jo 



2ci 



ln(cia;o + C2UJ0 + C3) 



C2 / 2ciu;o + C2- V^A 
\ 2ciu;o + C2 + 



(9.83) 



ln(c3) + 



where ci = m? and Cj = Ci{li,l2 



C2 



In 



C2 



C2 



2, 3. This expression, as it stands, is well 



defined and finite due to the presence of the test functions. To obtain a physically 
relevant result we have to remove this cutoff and therefore test the infrared behaviour 



of ( 9.83 ). We will show below that part of this expression becomes singular in the 
adiabatic limit. In the adiabatic limit in momentum space we let the test functions 
tend to a delta distribution in their argument. This will be done in the same way 
as in subsection 9.1, see ( |9.35 ). We scale the argument of the test functions with a 
parameter e in the integrals of ( 9.83 ) and consider the limit e — > 0. The functions 
C3 and A are homogeneous of degree 2 in the variables li and I2, i.e. 



C3(e/i,e/2,x) 
A{eli,el2,x) 



e^c^{li,l2,x) 
e'^A{h,l2,x) 



(9.84) 



while the function C2 is homogeneous of degree 1, i.e. C2{eli,el2,x) = ec2{h,l2, x). 
With this in mind let us now look at the various logarithms in ( 9.83[ ). The argument 
of the first logarithm contains the constant tti^Wq so it stays finite if e goes to zero. 
The constant in front of the second and the fourth logarithm has equal powers of 
e in the numerator as well as in the denominator so it tends to a constant. The 
argument of the second logarithm goes to one, so the logarithm itself goes to zero 
and the fourth goes to another constant. It remains to consider the third logarithm. 
This one becomes singular in the adiabatic limit. This can be seen as follows: 



da _ vr {2Tr)~^K^ 



d%d%g{h)g{l2rF2{+p-,)^ [ dx 

2ci Jo 



In 



{i0{pl)li - Pill cos {uj{pl)l2 - Pil2COS7)(l 



(9.85) 



+ finite terms. 
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If we concentrate on the singular part only we obtain by inserting the scaled test 
functions 



d^/i(i^/25(^i)5(^2)*ln|(tj(pj)/? - Pill cos + {uj{pl)l2 - Pil2COS7)(l - x)\ 

d'^hd'^l2ge{h)geih)*li'^\{^^{Pi)ll - PillCOS/?)x+ {uj{pi)l2 - Pil2COS7)(l - x)| 

= ^y" d^lid^l2go(^^go(^^^ ln|(u;(pi)/5' - pilicos/3)x+ - Pil2COS7) 
X (1 — x)| 

= j d'^kid'^k2go{ki)goik2)*^n\{uj{pi)eki - p\eki cos/?)x+ (uj{pi)ek2 - p\ek2 cosj) 

X (1 — x)| 
= (27r)4ln|e| +o(l) 

(9.86) 

where we've introduced ki = j. So we get a logarithmic divergence if e tends to 
zero. The origin of this singularity is the massless graviton propagator as can be 
clearly seen from the above calculation. Since all the other terms in ( |9.20 ) have the 



same structure of propagators we would get the same logarihmic divergence as in 
our example. Furthermore there can be no cancellation between these terms since 
they all have a positive sign. So we can omit the discussion of the other terms. 
This result is similar to the case of quantum electrodynamics |36] where a log- 



arithmic divergence in the differential cross section was obtained as well in the 
bremsstrahlung contribution to the scattering of an electron due to an external classi- 
cal source. As is well known from QED, the infrared divergence in the bremsstrahlung 
must cancel against contributions from radiative corrections, otherwise scattering 
theory would break down. 

9.3 Graviton and Matter Self-Energy 

In this section we discuss the infrared behaviour of radiative corrections to two 
particle scattering coming form the graviton self-energy and matter self-energy. Let 
us start with the graviton self-energy tensor. This is a fourth-order contribution to 
the S-matrix and the time ordered product T4(xi, . . . , X4) is given by 

Ti{xi, . . . ,X4) = in^ : ip{xi)^,^ip{xi),i3ip{x4)^,pip{xi)„ : 

-m^T^P" : ip{xi)^^ip{xi),fsip{x4)^ip{x4) : 

- m^jf^" : ip{xi)^ip{xi)ip{x4)lip{x^),^ : (9.87) 

-|- m^ry°^ry^°" : ip{xi)'^(f{xi)ip{x4)'^if{x4) : 

X D^ixi - a;2)6°^^^n(x2 - X3)xe^.ub'"'P'' {x^ - X4) 
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where the normahzation terms are already included in the graviton self-energy tensor 
n |26]. The S-matrix to fourth order is then given by 

54(5) = j d^xi . . . d^x^T^ixi, ... , X4)9ixi) . . . g(x4), g G cS(M^). (9.88) 

We consider the following initial and final states 

$i = / d'^pid^qiipi{pi)'>p2{qi)a{pi)^a{q{)'<n 

(9.89) 

^f= d^p2d^(B^{p2,q2)a{p2)^a{q2)^^. 

We are interested in the matrix element of S4^{g) between these states 

Sf, = {^f,SA{g)^i). (9.90) 

In the following we will consider the first term of T4 only, because all the other terms 
have the same propagator structure and differ from the first one in the powers of 
external momenta or powers of m only. We refer with an index (1) to the first term 
of (|9.87] ). Then the matrix element becomes 

Sf} = iK^^ j dVi d^qi d^P2 d^q2 d^xi . . . d^X4MPi)MQi)^if2, Q2) 

g{xi) . . . g{x^)D^{xi - :E2)5"^^^n(x2 - x^)x,^,h>^''"' {x^ - x^) (9-91) 
(r2,o(p"i)a(gi) : ^{xi)]^ip{xi),j3ip{xi)]pif{xA),a : a{pi)'^ a{q{)'^ 9) . 

We calculate the vacuum expectation value and write everything in momentum 
space. Then we obtain 

Sj^ = {27:)-Hk^ j dVi d'fi d^P2 dh~'2 d% . . . d%MP~i)MQi)^{P~'2, & 

g{h) . . . 5(/4)(l6a;(jri)^(Q-I)o.(p2)^(<r2))"'/'&"'''^6'^"^'^ 
Dq{12 + h + k-q2 + ^i)n(/3 + k-q2 + qi)\epiuD^{h - 92 + qi)P2Pl<l2'll 

+ D^ih + I3+U- q2+Pi)'n{h + k- q2+Pi)xe^uDo{U - 92 +Pl)P2 

+ {I2 + h + U-P2 + qi)'n{h + U-P2 + qi)xe^ub^ ik -P2 + qi)q2PiP2Qi 

+ b^ {h + h + lA-P2+ PlWz + k-P2+ Pl)\ef,ub^ ik -P2+ Pl)q2(llP2Pl 

5^^\h + ■.. + k- q2 + qi- P2+Pi)- 

(9.92) 

The adiabatic limit in the variable li is trivial, yielding just a factor (27r)^. Then 
we have to investigate the adiabatic limit in the remaining variables I2, ■ ■ ■ ,k- We 
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do this by introducing test functions with scaled arguments, see section 9.1. We 
consider the part depending on the variables h, ■ ■ ■ , ^4 in the first term of ( 9.92 ). 

d% . . . d%Uh) . . . ge{U)D^{h + k + h-q2+ qiW""' 
n(^3 + k-q2+ qi)xet.ub^''"'D^{k - 92 + qi)S^^^ {h + h + k - P2 + Pi - q2 + qi) 
£-12 j d^i^ . . . d'hUe-^h) ■ ■ ■ go{e-'h)D^{l2 + h + k-q2 + qiW^' 
n(^3 +k-q2+ qi)xef.ub^'''^D^{k - q2 + qi)S^^^ {k + h + k - P2 + Pi - q2 + qi) 
j d''k2 . . . d^hgoik2) . . . goih)D^{e{k2 + k3 + ki) -q2+ qi)b''''^' 

Il{e{k3 + k^) -q2 + qi)^,^h^'''^Do i^h - q2 + qi) 
5^"^^ {e{k2 + ks + k4) - P2 + Pi - q2 + qi) • 



(9.93) 



The graviton self energy tensor was explicitely calculated in [26| and it is given in 
momentum space by 



U{p) 



IT 



-d56pxPePfiPu - '^^^P^iPXPeVfiu + PtiPuriXe 



960(27r)5 

+ 1^2 p^{pxp^r]^y + pxPuVef, + PePfi-nXu + PePullXf,) 



- 162 p'^{rjXf,r]^u + VXyiletJi) + 118 p'^r]Xer]^u 
+ normalization. 



In 



-p^ - iO 
M2 



(9.94) 



where the normalization is a polynomial of degree four which has the same tensor 
structure as the one in (9.94). M > is an arbitrary mass scale. Now, according to 
( p. 93 ), the argument of the tensor 11 is given by 



P 



{e{k3 + k4) -q2+ qi)' 



2m^ 



2q2qi + 2eik3 + A:4)(-g2 + ^i) + e^iks + k^) 



(9.95) 



since the momenta qi and q2 have to be on the mass shell. We see that in the adiabatic 
limit e ^ this goes to 2(m2 — q2qi)- The Feynman propagator in momentum space 
is given by 



bE{p) = -{2^)- 



1 



p^ + iO 



(9.96) 



Again we observe that in the adiabatic limit the arguments of the two propagators 
in ( p. 93 ) approach 2{im? — q2qi)- Therefore, ignoring the tensor structure of 11, we 
can argue that in the adiabatic limit we obtain 



/ 



d k2 ■ ■ - d A;4 . . . = In 



-2(m2 - g2gi) 
M2 



+ o{l) 



(9.97) 
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where we have used that 

^ d'^kgoik) = {27Tf. (9.98) 



This result is finite due to the mass m in the numerator of the logarithm. So we 
conclude, that the adiabatic limit for radiation corrections coming from the graviton 
self-energy tensor is infrared finite. 

Now we come to the matter self-energy. As in the previous case we restrict ourself 
to a typical term of the time-ordered product which shows all essential features of 
the infrared behaviour. Furhtermore it is enough to consider the matter self-energy 
function in a scattering process where one massive scalar particle is scattered at a 
fixed external gravitational field. The time-ordered product to third order has the 
form 

T3(xi,X2,X3) ~ iK^ : ip{xi)^ip{x3) : E(xi - X2)D^{x2 - X3)/i^^'*(x3) (9.99) 

where h^^^ is a classical external gravitational field and S is the matter self-energy 
function. The S-matrix to third order is then given by 



S3{g)= J d''xi...d''x3T3{xi,X2,xs)g{xi)g{x2)g{x3), <? G 5(M*). (9.100) 

Then we want to calculate the S'-matrix element 

Sfi = {^f,Ss{g)<^i) (9.101) 

where $j and are one-particle wave packets in momentum space. The S-matrix 
element in momentum space then takes the form 

Sr^=-^ [ -^P=-^3=';p{p)i!{q) [ d%...d%g{h)...g{k) 

±{-h + q)D^{-h -h+ gyh'^^'i-h -h-h + q-p). 

In this expression we can perform the adiabatic limit in the variable ^3, since it 
appears in the argument of the external source only. We set this variable equal to 
zero in /i^^* and obtain a factor (27r)^ from the integration. To study the singularity 
structure we look at the integrals w.r.t. h and h- 

J d\d% ...= I d^kid%go{ki)go{k2)±{-eki + q) ^^^^ 

D^{-eki - ek2 + q)}f''\-eki - ek2 + q - p). 



The matter self-energy function S including it's normalization was calculated in |24], 
so we quote here the result in momentum space only. 




iird^p^ — m 



2\ 



(9.104) 
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where C2 is a normalization constant and T is given by 



-m?TT 



2(27r) 



4 ■ 



(9.105) 



Taking only the leading singularities in the Feynman propagator D^{q — e{ki + ^2)) 
it becomes 

m,-e(h+k,)) = i2,r'^-^^-l^-^^. (9.106) 

The self energy function has singularity structure of the following form 

2eqki 



-2eqki In 



+ 0(1). (9.107) 



Then we obtain an overall logarithmic divergence 

/ d%d^k2 . . . = In |e| / d^ki<fk2m{ki)m{k2) ^ „ , + o{l). (9.108) 

J J 2q{ki + k2) - iO 

since the factor e infront of the logarithm is canceled by the Feynman propagator. 

In the theory of massive scalar matter coupled to gravity there remains to discuss 
the so called vertex-function which is again a third-order process. Due to the self- 
coupling of the gravitational field, which was extensively studied in [^ ] there exists 
vertex functions of different type. This is not the case in QED. Since all infrared 
divergences calculated so far have the same form ~ Ine, it is quite plausible that 
they exactly cancel as they do in QED. In this case Weinberg's short argument 
would be right. For the full proof one needs the various vertex corrections. This is 
beyond the scope of the present work. 



„Alle Anfange sind dunkel. Gerade dem Mathematiker, der in seiner 
ausgebildeten Wissenschaft in strenger und formaler Weise mit seinen 
Begriffen operiert, tut es not, von Zeit zu Zeit daran erinnert zu werden, 
dafi die Urspriinge in dunklere Tiefen zuriickweisen, als er mit seinen 
Methoden zu erfassen vermag. Jenseits alles Einzelwissens bleibt die Auf- 
gabe, zu begreifen. Trotz des entmutigenden Hin- und Herschwankens der 
Philosophic von System zu System konnen wir nicht darauf verzichten, 
wenn sich nicht Erkenntnis in ein sinnloses Chaos verwandeln soil." 

Hermann Weyl, 1918 in Raum ■ Zeit ■ Materie 
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Appendices 



In the subsequent appendices A and B we give the exphcit divergence forms for the 
various types in dqTi. 

A Divergences for Types ^4, 5, C, H and K 



Here we give the unique divergence expressions for dqTi. 
1. Type A: We calculate the expression d^Ty^ explicitly: 

a^ff/'f = di u'^,^h'ShP'' + d2 u'i,hP;,^hf^ + ds U^^hP^hP^ 

The constants di, . . . , are given by 

di:=ci + C5, d2:=ci+C7, := ci + ce + cy, 
di ■■= C2 + C5, ds := 2 C2 + cg 

From first order gauge invariance we obtain 



(A.l) 



(A.2) 



1 

di = 0, d2 = 0, ds = —Ms, = 2 ^8' ^5 = 



The coefficient matrix Ma of ( |A.2| ) is in GL(5,Z). We invert these equations and 
obtain 

ci = di + ^{d2-d3 + d5-2di) (A.4) 

C2 = ^ (4 - + ^^2) (A.5) 

C5= ^(f^3-d2-4 + 2(i4) (A.6) 

C6 = ds- d2 (A. 7) 

C7 = d4 - + - {d2 +d3- ds) (A.8) 

These equations give, together with (|A.3|) , the desired divergence for dgTilxypeA- 

p,r 
1/1 



2. Type B: We calculate the expression d^T'^j^ explicitly: 



^M^i/i = ^6 ut^,Ji^h + dj u^^h,^,^h + ds u^h^^h^ 

The constants ) • • • > ^^lo are given by 

d&:=h + ci2, d7:=C8 + ci4, := cs + C13 + cm, 
d9:=C9 + Ci2, dio:=2c9 + ci3 



(A.^ 



(A.IO) 
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From first order gauge invariance we obtain 
de = 0, dj = 0, ds '- 



ia2, dg = -a2, dio = 



(A.ll) 



The coefficient matrix Mb of ( [A.IO ) is in G-L(5,Z). We invert these equations and 
obtain 



cs = + ^ {dio -2dg-ds + dj) 



eg 



(dio - 4 + d^) 



Cl2 = - (4 + 2 dg - dio - dj) 
ci3 = ds-dr 

ci4 = dg - (ie - ^ (c^io - ds - dj) 



(A.12) 

(A.13) 

(A.14) 
(A.15) 
(A.16) 



These equations give, together with ( A.ll| ), the desired divergence for dqTilTypeB- 
3. Type C: We calculate the expression dyT'^j-^ explicitly: 



a^ff;f = dii uZ,^hZ^'h + di2 ulhZ^t^h + di3 ulhZ^'h^ 

+ di4 uZ,f,h''f'h,^ + di5 uZhZ^K + dw uZh'^f'h,^,, 
+ dir uZhZf^h,, + di8 n"/i^^/i,, + dig 

The constants dn, . . . , dig are given by 



(A.17) 



dii 
di4 
di7 



C15 + C21, di2 := Ci5 + C24, di3 := C15 + C25 + C26, 
CI6 + C22, dl5 := Cl6 + C23 + C24, die := CI6 + C25 
C17 + C21+C22, di8 := Ci7 + C23, dig:=Ci7 + C26 



(A.18) 



From first order gauge invariance we get 
dii = 
rfi5 = 



iae, di2 = 0, di3 
i 



■as, d 



14 



■ as, 



2 2 

04, die = 0, di7 = 0, di8 = 0, dig = 



(A.19) 



The coefficient matrix Mq of ( |A.18 ) is in GL(9, Z). We invert these equations and 
obtain 



1 



Cl5 = 7^ (^11 + dl2 + di4 - di5 - di7 + dis) 



ci6 = 7^ {dn - di3 + di4 + die - dn + dig) 



cn = 77 {di2 - di3 - di5 + die + dig + dig) 



(A.20) 
(A.21) 
(A.22) 
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C21 = 


i 

2 


(c^ii - 


di2 - du + di5 + di7 - 


dis) 


(A.23) 


C22 = 


1 

2 


{-du 


+ dn + (ii4 - diQ + di7 


- dig) 


(A.24) 


C23 = 


1 

2 


{-di2 


+ c^i3 + di5 - die + dis 


- dig) 


(A.25) 


C24 = 


1 

2 


{-du 


+ di2 - du + di5 + (ii7 


- dis) 


(A.26) 


C25 = 


1 

2 


{-du 


+ di3 - du + die + di7 


- dig) 


(A.27) 


C26 = 


1 

2 


{-di2 


+ dis + (ii5 — die — di8 


+ dig) 


(A.28) 



These equations give, together with ( A.19| ), the desired divergence for dqTilxypeC- 



4. Type D: We calculate the expression d^T'^'P explicitly: 



a^Tfyf = d2o + d2i ulh%%h''^ + d22 uZhl^h^^ 

+ d23 U^h'^-hl^ + d24 + (i25 uZh'^'^hl^, 

+ d2e u^hZ^hl^ + d27 u^K^,,hl^ + d28 
The constants ^20, • • • , (^28 ^ire given by 

^20 := C27 + C33, ^21 := C27 + C36, ^22 := C27 + C37 + C38 
^23 := C28 + C34, ^24 := C28 + C35 + C36, ^25 := C28 + C37, 
C^26 := C29 + C33 + C34, ^27 := C29 + C35, t^28 := C29 + C38 

From first order gauge invariance we obtain 



(A.29) 



(A.30) 



d20 
d24 



- o,g, d2i — 0, d22 



"2 '^10; d23 



ag, d: 



25 



0, d 



26 



0, d 



'27 



-iai2 
0, d28 = 



(A.31) 



The coefficient matrix Md of ( A.SOj ) is in GL(9,Z). We invert these equations and 
obtain 



C27 

C28 

C29 
C33 
C34 
C35 



{d20 + d2l + d23 - ^24 - d26 + ^27) 
{d20 - d22 + d23 + d25 " C?26 + C^28) 
(^21 - d22 - ^24 + d25 + (^27 + (^28) 
{d20 - d2l - d23 + ^24 + d26 - ^27) 
(-^20 + d22 + d23 - d25 + ^26 " C?28) 

6^28) 



2 (-'^21 + C^22 + C?24 - (^25 + <^27 



(A.32) 
(A.33) 
(A.34) 
(A.35) 
(A.36) 
(A.37) 
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21 



C36 = ^ (-^20 + d 
C37 = ^ (-^20 + f^22 



C?23 + d2A + ^26 - ^^2?) 
C?23 + f^25 + ^26 - C^28) 



1 



C38 



-d2l + d22 + d2A - d25 - ^27 + C?28) 



(A.38) 

(A.39) 
(A.40) 



These equations give, together with ( A.31| ), the desired divergence for dqTilTypeD- 



5. Type H: We calculate the expression d^Ty-^ explicitly: 



dfj.T^l" = dso ut^^,^u1u" + dsi u^u"^,^u°' + ^82 uf^uZu'^^ 



+ ds3 u1,^u°'uZ + 44 u'liL^^uZ + ds5 

+ ds6 u^^uZuZ + ds7 u^uZ,^uZ + (i88 u^uZu"^,^ 



(A.41) 



The constants dgo, . . . , ^88 are given by 

C84 + C87, dsi := C84 + C90, ds2 ■= C84 + C91 + C92, 



dso 

d83 

dse 



C85 + C88, 44 := C85 + C89 + C90, ^5 := C85 + C91 , 
C86 + C87 + C88 5 ^7 := C86 + C89, 48 := C86 + C92 

From first order gauge invariance we obtain 



(A.42) 



d 



80 



=0, 41 = 0, di 



82 



519 5 



d: 



83 



0, 4 



'19 



(A.43) 



45 = - 2^2, 46 = "2 (^4 + &19), 47 = 0, 48 = --63 

The coefficient matrix M// of (|A.42| ) is in GL(9, Z). We invert these equations and 
obtain 



C84 = 


1 

2 


^4o + 4i + 


43 — 44 — 46 + 47) 


(A.44) 


C85 = 


1 

2 


^4o — 42 + 


43 + 45 — 46 + 48) 


(A.45) 


C86 = 


1 

2 


^4i — 42 — 


44 + 45 + 47 + 48) 


(A.46) 


C87 = 


1 

2 


^4o — 4i — 


43 + 44 + 46 — 47) 


(A.47) 


C88 = 


1 

2 


^— 4o + 42 


+ 43 — 45 + 46 — 48) 


(A.48) 


C89 = 


1 

2 


^— 4i + 42 


+ 44 — 45 + 47 — 48) 


(A.49) 


C90 = 


1 

2 


^— 4o + 4i 


— 43 + 44 + 46 — 47) 


(A.50) 


C9I = 


1 

2 


4o + 42 


— 43 + 45 + 46 — 48) 


(A.51) 


C92 = 


1 

2 


^— 4i + 42 


+ 44 — 45 ~ 47 + 48) 


(A.52) 
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These equations give, together with ( A.43| ), the desired divergence for dgTilTypeH- 

/I 



6. Type K: We calculate the expression dpif'^',^ explicitly: 



dfiT^/f = dio2 u'^^^^ulu'' + (iio3 ulu'l,^u'' + dio4 u'^^u^u'^^ + dios li^.^u^u^ (A. 53) 

The constants dio2, • • • , di05 are given by 

diQ2 '■= cio8 + Clio, t^ios := cio8 + C112, dio4 := — cios + cuo + cm, dies := cm 

(A.54) 

From first order gauge invariance we obtain 

dio2 = 2 ^18, dio3 = 0, dio4 = "2 ^13, di05 = 2 ^17 (A.55) 

The coefficient matrix Mk of ( |A.54| ) is in GL(4, Z). We invert these equations and 
obtain 

Clio = ^ (rfi02 - c?i04 + c^ios) (A.56) 

C112 = ^ (dio2 + fiio4 - f^ios) (A.57) 
C113 = c^i05 (A. 58) 

C114 = dios + ^ {-dio2 + dio4 - dies) (A.59) 



These equations give, together with ( A.55| ), the desired divergence for dqTilTypeK ■ 



B Divergences for Types E, F,G and J 

Here we calculate the explicit divergence forms in terms of the coupling parameters 
, 012, bi, . . . ,^21 for the types E, F, G and J. In contrast to the other types 
the system of equations between the Cj and dj are no longer invertible in a unique 
way. There are some ambiguities, if we express the Cj in terms of the dj. Let us 
begin with type E. 

1. Type E: Let Me G Mat{l8 x 15, Z) the coefficient matrix of the system (JT^. 
We have to consider the equation 

Me ■c^ = d^ (B.l) 

The general solution of this equation is the sum of an arbitrary solution and the 
general solution of the corresponding homogeneous equation 

Me • = (B.2) 
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The matrix Me has rank{ME) = 12. So there are three free parameters Ai, A2, A3 E 
C in the solution of (|B.2|). We obtain 



Cq (Ai, A2, A3) — — A2 + A3, — Ai + A2 — A3, Ai — A2 — A3, — Ai — A2 + A3, 

— Ai — A2 + A3, — Ai + A2 — A3, — Ai + A2 — A3, Ai — A2 — A3, 
Ai — A2 — A3, Ai + A2 — A3, Ai — A2 + A3, — Ai + A2 + A3, 2 Ai, 

2A2, 2 A3) 

(B.3) 

A special solution of equation ([B.l|) is given by 



2 (^37 + dm — dzo) , ^29 — C?31 — (^44 + (^42 — 2 (^39 — d^j — d^o) , 
1 1 

^32 + dse — ^33 — d44 — - [dsg — ^37 — d^o) , ^29 — d^ — - {dsg — ^37 — ^40) , 

1 1 

^32 — ^44 — - (c^39 — '^37 — d^o) , —^32 + C^33 + ^44 + - (c^39 — d^l — C^4o) , 
^ (d39 + d^j - (i4o) , ^ (^40 + (i39 - ^37) , ^31 - ^29 +^44 + ^ (c?39 " C^37 

— d^o) , ^44 + - {d^g — d^i — dio) , ^44 + dsi — ^29 — ^42 + ^45+^ ('^39 

- fi37 - d4o) , -c^32 - dsQ + ^33 + d^Q + du + - [dsQ - d37 - d^o) , 0, 0, 



The general solution of (B^) is then given by 



(B.4) 



(B.5) 



With the equations (|7.4| ) we can write the expression dqTilTypeE as a divergence. 
The two parts of the solution to (B.l), cf and c^, correspond to a Q-vertex of the 
form 



^1/1 ~ ^1/1 + ^1/1 



/I 



1/1 



(B.6) 



where the distinguished part T^j^ is given by (6^) with coefficients (B^). They are 
uniquely determined by the parameters of the theory, see ( |7.4| ). The homogeneous 
part Ty^'^ is given by ( |6.5D with coefficients (|B.3| ) and one observes that it can be 
written in the form 



rf;f'°= (Ai + A2-A3)aprr' 



pp,E 



+ (Ai-A2 + A3)apr^ 
+ (Ai-A2-A3)apr3^''' 



(B.7) 
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with 



j,^P,E ^ ^^f^»j^f^P^ _ u'^hlPh^" (B.8) 
they are antisymmetric in their indices, i.e. 

rpPP,E ^ _ri.pp,E 1^2,3 (B.9) 

The homogeneous part Ty^'^ can be written as the divergence of an antisymmetric 
tensor which is independent of the parameters of the theory. 

2. Type F: In analogy to type E we first calculate the expression dfjTy^ . This 
expression is of the form 

61 

^/^^V'f = E '^^ d^dpd„\uVh'''' (B.IO) 

i=47 

Here the three derivatives are distributed among fields in all possible combinations. 
The new constants di are defined as follows 



^47 

dm 



C57 + C58 + 2 C63, ^48 := C56 + C57 + Cq-j, d^g := C57 + C59 + C66, 
C58 + C62 + C67, ^^51 := C55 + C58 + C66, ^52 := C59 + C62 + C64, 

C55 + C59 + 2 C65, := C54 + C66 + Ceo, t^ss := ceo + cei + C64, (B.ll) 

Ceo + C65, C?57 := C54 + C61 + C67, ^58 := C61 + C68, 

C56 + C62 + 2 C68, f^60 := C54 + C63, C^61 := C55 + C64 + C56 



From first order gauge invariance we obtain 

d47 = - i(^^ag + ai2^, di^ = -ihi, ^49 = -^09, 
^50 = - ■if^ aio + &2y 41 = 0, d52 = -ib3, 43 = 0, 



.2 

c?54 = 1(^05 + - ag + ai2^, d55 = -ib4, 46 = ^(206 + 09 + 012), (^•^^) 

4? = 2 (2 "4 + 09 + oio) , '^ss = 2 (2 ^3 + oio + an - 2 65) , 
49 = - ibe, 4o = 0, 4i = 

Let Mp e Mat{15 x 15, Z) the coefficient matrix of ( p?Tl| ). Then we determine the 
general solution of 

Mf-c^ = (B.13) 
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where G C"^^ and d-^ G C^^ are the column vectors with components (C54, . . . , cgs) 
and ((i47, . . . ,^61) respectively. The matrix Mp has rank{Mp) = 11. The general 
solution of the corresponding homogeneous system 

Mp-c^ = Q (B.14) 

is labeled by 4 independent parameters Ai, . . . , A4 G C and is given by 

Cq{\i,\2, A3, A4) = ^— A2 + A3, Ai + A3 — A4, — Ai — 2 A3, Ai — A2 + 2 A3, 

— Ai — A2, — Ai — A3 — A4, — A4, —A3, Ai, A2 — A3, A3 + A4, 
A4, A2 — A3 + A4, A2, A3 



(B.15) 



A special solution to ( |B.13| ) is given by 



C^56 + C^57 — C^55 — 2 ('^^3 + C?59 — d52 — C^6l) 1 ^ ('^53 — d^g — ^52 + del) , 

(^59 1 dig + d^j — (iss — d^i — d^s — d^g + (igi + 2 d^Q, ^47 — d^g + ^52 

+ d^j — (iss — 2 deo + ^54, ^ ('^53 + ^59 + ^52 — (igi) , ^56, —^56 + ^55 
+ 2 (^53 + ^59 — (^52 — del) , 0, —^56 — ^57 + ^55 — deO + ^ ('^53 + dsg 

— d52 — del) , — - (d53 + dsg — d52 — del) , 0, d54 — 2 dse — d^j + dss 



+ 2 (^53 + d59 - d52 - del) , 0, 



(B.16) 



The general solution to (B.13) is then given by 



F , F 



(B.17) 



With the equations ( B.12 ) we can write the expression dqTilxypeF as a divergence. 
According to the two parts of the solution to ( [B.15 ), and Cq" we can represent 
the Q-vertex as a sum 



-'l/l ~ -'l/l +-'1/1 



(B.18) 



The distinguished part Ty^'^ is given by ( |6l^ ) with coefficients ( B.16 ). They are 
uniquely determined by the parameters of the theory, see (B.12). The homogeneous 
part Ty^'*^ is given by ( |6.6| ) with coefficients ( p.l5| ) and it has the form 



(B.19) 
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with 







uPhp^Kp^ 


rj,jJ.p,F _ 


uPhP^h"" - 


uPhP^h"" 


rpfJ.p,F _ 

^3 — 


uP,h'"'h''P - 


uP^hP'^h'"' 


rplJ.p,F _ 


uPhP'^h':^ - 


uPhP'^h"^ 



(B.20) 



They are antisymmetric in their indices, i.e 

rj.^.a,F ^ _j.a,,,F^ i = 1, . . . , 4 (B.21) 

The homogeneous part Ty^''^ can be written as the divergence of an antisymmetric 
tensor which is independent of the parameters of the theory. 
3. Type G: We calculate d^TP'^^ . This has the form 



d.TPj^ 



79 

E 

i=62 



didadpdAuPhP^h 



The new constants di are defined by 



(B.22) 



d&2 
d^b 
des 
dn 
dji 
dn 



C72 + C78, 

C73 + C78, 

C74 + C79, 
C75 + C79, 

C76 + C80, 
C7I + C80, 



d&z 
dm 
dm 

d72 
d75 

djs 



C70 + C72 + C81, 
C73 + C74 + C81, 

C7I + C74 + C83, 
C75 + C76 + C83, 
C70 + C77 + C83, 
C7O + C79, (^79 



^64 
(^67 
^70 
^73 
^76 



C72 + C77 + C82, 
C7I + C73 + C82, 
C69 + C75 + C81, 
C69 + C76 + C82, 
C77 + C80, 



(B.23) 



C69 + C78 



From first order gauge invariance we obtain 



^62 

dee 
djo 

d72 

diA 



iae, des = -ibr, 
1(^^0.4 + bs^, der 



i( ai + 06 + - 07 



d 
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^68 = ■ 

-ibio, 



^65 — — 
-ibg, dm 



0, 



i /I 
i{ai - bii) + - (04 + 05 + 07), (i73 = i(2a2 + - as + 



0, d 



75 



2 

-^^12, 



76 



0, d-r-j = 0, 



78 



2 

0, (^79 







(B.24) 



Let Mg e Mat{18 x 15, Z) be the coefficient matrix of ( |B.230 . We determine the 
general solution of 



Mg • c' 



G 



(B.25) 
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where c*^ G C"^^ and d*^ G C"^^ are the column vectors with components (cgg, . . . , cga) 
and (^62, • • • ,^79) respectively. The matrix Mq has rank{MG) = 12. The general 
solution of the corresponding homogeneous system 

Mg-c^ = Q (B.26) 

is labeled by three independent parameters Ai, A2, A3 € C and is given by 

Co'(^i) ^2, A3) = — Ai — A2, A2 — Ai — A3, Ai — A2 — A3, A3 — Ai — A2, A3 

— Ai — A2, A2 — Ai — A3, A2 — Ai — A3, Ai — A2 — A3, Ai — A2 

— A3, Ai + A2 — A3, Ai — A2 + A3, A2 + A3 — Ai, 2 Ai, 2 A2, 2 A3^ 

(B.27) 

A special solution to ( p.25| ) is given by 

= (^73 - + ^70) , ^ (dee + ^67 - c^eg) - ^65 - ^^64 + ^62 + ^75, 

1 11 

- (^67 — (ige + (ieg) , '^62 — d^^ — - (deg — c^ee — ^^67) 1 2 (^66 + d^i — ^eg) , 

^ (deg + dQQ - der) , ^ (^72 - '^73 + djo) , ^ (^72 + ^73 - ^70) , ^5 + c^64 

- dG2 + - {dm — dae — d^-j) , dgs + - [dag — dee — ^67) , (^65 + '^64 — ^62 

- d75 + drs + ^ {dm - d^-j - dee) ,^74 + ^ (dee - c?67 - deg) , 0, 0, 0^ 

(B.28) 



The general solution to (B.25) then reads 



-G 



(B.29) 



With the equations (B.24) we can write the expression dQTi\TypeG as a divergence. 
According to the two parts of the solution to ( B.25 ), and Cq', we can represent 
the Q-vertex as a sum 



^1/1 ~ ^1/1 + ^1/1 



/I 



1/1 



(B.30) 



where the distinguished part T^'^ is given by (6/7) with coefficients ( B.28| ). They 
are uniquely determined by the parameters of the theory, see ( B.24[) .The homoge- 
neous part T^j'^'^ is given by ( |6.7|) with coefficients ( B.27 ) and one observes that it 
can be written in the form 



^1/1 



(Ai + A2-A3)aprr'^ 
+ (Ai-A2 + A3)apr2^'''^ 
+ (Ai-A2-A3)apr3'^^'^ 



(B.31) 
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with 



rppp 



G 



(B.32) 



rppp 



G 



uPh^'^h^ - u^hP^'h,, 



They are antisymmetric in their indices, i.e. 

rpp.p,G rpPP,G 



1,2,3 



(B.33) 



The homogeneous part Ty'^'^ can be written as a divergence of an antisymmetric 

tensor which is independent of the parameters of the theory. 

4. Type J: We calculate the expression dfjT'^j'^. This has the form 



101 



dpTPjl = Y,did^da,dp\uP 



u'^uP 



(B.34) 



The new constants di are defined by 



dm 



— C96 + Cioi — Ci02, dgo := — Cge + Ci04 — Cioo + Cl08, 

C96 — Ci03 + Ci07, ^92 = Ci03 — Cioi, C?93 := " C97 + C107 — C99, 

C97 — Cioo + Ci05 + C108, dg^ := C97 — C102 + C106, 

Cl07 — C102, dgj := Ci03 — C98 — C106, ^98 := C98 — Ci04 + Ci05, 

C106 — C99, dioo := C98 + C99 — Cloi, dioi := cioo + cio8 



From first order gauge invariance we obtain 
^89 = 



- {b3 + bis) , dgo = 0, dgi = - - {bi - bis) , 



dg2 
dm 

0?98 



2 (^2 + bn) , ^^93 = 2 ^18 + ^^12' dg/i 



0, 



- (61 + 64 + 2 67 + 613) , dgQ = 0, dgj = - 62 + i68, 

i 



■bs-ibg, di 



'99 



bn + ib 



ii> 



100 



0, d 



101 







(B.35) 



(B.36) 



Let Mj G Maf (13 x 12, Z) be the coefficient matrix of (|B.35D . Then we determine 
the general solution of 



Mj-c 



(B.37) 



where G C"*^^ and d"^ G C^'^ are the column vectors with components (C96, . . . , C107) 
and ((^89,... ,(iioi) respectively. The matrix Mj has rank(Mj) = 9. The general 
solution of the corresponding homogeneous system 



(B.38) 
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is labeled by three independent parameters Ai, . . . , A3 G C and is given by 
c^(Ai,... ,A3) 



Ai, A3 — A2, Ai — A2 + A3, A2, 0, 
Ai + A3, A3, Ai + A3, Ai, A2 — A3, 
A2, A3 



(B.39) 



A special solution to ( |B.37|) is given by 

""^ — ' dioi — dgo — 2 {dgs — dlOO — dgg — dsg), —dgi — dgj — dioo — dgg — dgg + ^95^ 
dlOl — dgo — dgi — dgj, —2 dgs + 2 dloo + dgg, dioi , 

C^lOl — dgo — dgi — dgj — 2 dgs + diQQ + dgg, —dgi — dgj — diQQ — dgg — dsg, 
dlOl — dgo — dgi, 0, dioi + dgo + dgi + dg-j + dgs, 0, 



The general solution to ( B.37 ) is then given by 



(B.40) 
(B.41) 



With the equations ( [B.36| ) we can write the expression dQTi\Typej as a divergence. 
According to the two parts of the solution to ( B.37 ), and Cq, we can represent 
the Q-vertex as a sum 



^1/1 
,J,S 



^1/1 



^1/1 



(B.42) 



where the distinguished part T^J-^'"^ is given by ( |6.9| ) with coefficients (B.40). They 
are uniquely determined by the parameters of the theory, see ( p. 36 ). The homoge- 
neous part Tl^j'^'^ is given by ( |6.9D with coefficients (|B.39D and one observes that it 
can be written in the form 



rf,'/'0 = (A2 - x^)dpT^p^' + Ai^^rr-" + \zdpT^P'- 



-1/1 



(B.43) 



with 



rpfJ.p,J 










- uPpU'^uP 


rpllp,J 




- uii^u^'uP 



(B.44) 



They are antisymmetric in their indices, i.e. 

TPP'-^ = -TPP"^, i = 1, 2, 3 (B.45) 

The homogeneous part T^'^''^ can be written as a divergence of an antisymmetric 
tensor which is independent of the parameters of the theory. 
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